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ABSTRACT 


The orientation dependence of the de Haas—van Alphen periods in tin 
single crystals has been studied in high pulsed magnetic fields. As many as 
seven periods have been found at some orientations, indicating that the 
Fermi surface is much more complicated than the 3-ellipsoid model suggested 
by Fawcett. The general degree of complexity of the results is much the 
same as that predicted by the free-electron model, but serious discrepancies 
are found when a detailed comparison is made. Effective mass values have 
been determined at the principal orientations from the temperature depen- 
dence of the amplitudes; these are compared with the results of cyclotron 
resonance experiments and with the free-electron orbit masses. 


§ 1. InTRODUCTION 


Tue period, P, of the de Haas—van Alphen oscillations in a metal single 
crystal at low temperatures is related to the extreme (maximum or 
minimum) area of cross section, ,, of the Fermi surface normal to the 
magnetic field by the Onsager—Lifshitz and Kosevich relation : 


PED nolchil ae 8 at le GAL) 


(., in k-space), and from the temperature dependence of the oscillation 
amplitudes one can obtain the corresponding cyclotron (orbit) effective 
mass, m*, which is defined by: 
Mee ca (OM /OL oN. wer ta ood ay 4a (2) 
In principle a study of the variation of P and m* with orientation is 
sufficient to determine the shape of the Fermi surface and the Fermi 
velocity at all points on the surface, and for a discussion of the present 
theoretical and experimental situation, the reader is referred to the 
review articles by Shoenberg (1957, 1960). 
The orientation-dependence of the de Haas—van Alphen effect in metallic 
tin has been studied at magnetic fields up to 15 ke principally by Shoenberg 
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(1952) and Berlincourt (1952), whose findings were, for the most part, 
consistent with a Fermi surface consisting of two small ellipsoids oriented 
along the tetragonal axis and capable of accounting for a total of about 
10-2 carrier per atom at most. Subsequent preliminary studies using 
impulsive fields of up to 100ke (Shoenberg 1953) revealed a rich variety 
of short-period de Haas—-van Alphen oscillations associated with the 
major portions of the Fermi surface. A detailed orientation-dependence 
study of the periods in lead using the impulsive-field technique has been 
successful in determining the Fermi surface in this metal (Gold 1958; 
hereinafter referred to as (1)) and in this paper we describe such a study 
of metallic tin. 

The high-field de Haas-van Alphen results turn out to be very com- 
plicated and are, in fact, of the same order of complexity as predicted by 
the simple free-electron model. We shall make use of this simple model, 
which proved helpful in the interpretation of the lead results, in an attempt 
to relate the observed results to the Brillouin zone scheme. However, 
we shall conclude that further information, in particular from de 
Haas—van Alphen studies at still higher fields and from ultrasonic attenua- 
tion, cyclotron resonance and magneto-resistance experiments, is 
necessary before all of the geometrical features of the Fermi surface can be 
established with reasonable certainty. We have also determined the 
effective masses corresponding to principal orientations; these are com- 
pared with the results of cyclotron resonance experiments and with the 
orbit masses predicted by the free-electron model. 


§ 2. EXPERIMENTAL RESULTS 


The experimental method is the same as that described in (I) and in 
the review article by Shoenberg (1957), so that no further details will 
be given here. Since the high-field results are so complex, with up to 
seven periodic terms superimposed on one another at some orientations, 
it was essential to use the resonance method in order to sort out the 
various terms; the presence of a background of the other terms 
nevertheless resulted in a greater uncertainty of period measurement 
than was the case for the simpler lead results (I). The measurements 
were carried out on single crystal wires of length 8mm and diameter 
0-55mm; these were prepared from Johnson Matthey spectroscopically 
pure tin and grown with a principal crystallographic axis along the wire. 
As was explained in (I), the absolute orientation of the crystal axes 
relative to the magnetic field was obtained finally from the symmetry 
of the period-orientation results. The period measurements were carried 
out at 1-2°K and using field pulses rising up to 80ke in about llms. A 
few measurements were also made along [001] and [110] in a field of 
150ka@; this became possible (using a specially designed field coil) only 
after the bulk of the measurements had been completed. The higher 


field greatly increased the observed amplitudes but revealed no new very 
short periods. 
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Variation of period with orientation of magnetic field in a (100) plane. Upper 
half: long-period oscillations; lower half: short-period oscillations. 
©, A, specimens with [001] axis along wire; @, specimen with [100] 
axis along wire; +, specimen with axis midway between [100] and 
[001] along wire. — - — -, Shoenberg (1952); ——-, Berlincourt (1952). 
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Fig. 2 
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Variation of period with orientation of magnetic field in a (110) plane. Upper 
half: long-period oscillations; lower half: short-period oscillations. 
A, specimen with [001] axis along wire: Cj, specimen with [110] axis 
along wire. —-— ., Shoenberg (1952). i 
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Figures 1, 2 and 3 show the variation of the periods with orientation of 
magnetic field in the (100), (110) and (001) planes, respectively. The 
results fall, broadly speaking, into two groups: long-period oscillations 
with P~10-7@- (upper diagrams) and short-period oscillations with 
P~10-%a@™ (lower diagrams). Also shown in the figures are smooth 
curves drawn through the long-period data obtained earlier by Shoenberg 


Fig. 3 
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iati f period with orientations of magnetic field in the (001) plane. 
ig pee half : long-period oscillations; lower half: short-period oscilla- 
tions. @, specimen with [100] axis along wire; [J], specimen with 

[110] axis along wire. —-—., Shoenberg (1952). 
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(1952) and Berlincourt (1952) using the couple method at fields up to 
15ke¢, and it is satisfactory that the corresponding high-field results 
agree with the earlier low-field data. In table 1 we give values of the 
period and of the associated orbit effective mass at the principal orien- 
tations [001], [100] and [110]. 


Fig. 4 
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Because of the considerable scatter of the experimental points and the 
rapid orientation dependence of amplitude (see fig. 4), it is not possible 
to decide unambiguously how to join up the experimental points into 
curves, without appeal to some theoretical model. For instance it is 
not clear whether the curve L in fig. 1 is a separate branch existing over 
only a limited angular range or whether it is an extension of E or of F, 
and there are clearly many similar ambiguities. For lead (I) certain 
features of the Fermi surface could be inferred by straightforward inspec- 
tion of the data and the surfaces were then fitted into the Brillouin zone 
scheme using the free-electron model as a guide; however, we now find 
that for tin some, theoretical model is necessary to suggest the very 
nature of the surfaces themselves. 


§ 3. THe Free-ELectron MopEL 


The simple free-electron model has been remarkably successful in 
accounting for the geometrical features of the Fermi surfaces in the 
polyvalent metals lead (I) and aluminium (Heine 1957, Harrison 1959, 
1960) and it now appears to have a good theoretical foundation (Harrison 
1960, Cohen and Heine 1960). We now appeal to this model for help 
in understanding the complex period variations in figs. 1, 2 and 3. 

The primitive Brillouin zone for the white-tin lattice is shown in fig. 57. 
Since the primitive zone always accommodates two electrons per smallest 
unit cell of the real lattice and since for tin the smallest unit cell contains 
two atoms, it follows that the zone can accommodate one electron per 
atom. Assuming tin to have four valence electrons per atom, these can 
be accommodated within a sphere in k-space of volume 8(27/a)%(a/c) em=3 
and hence radius 1:518 (27/a)cm—!, where a=5-82A and c=3-18A are 
the lattice constants. This sphere completely encloses the primitive 
(first) zone so that the Fermi surface must lie entirely in higher zones in 
the extended zone scheme. In order to make comparison with experi- 
ment it is necessary to re-map the various portions of the surface in the 
higher zones into the primitive zone by translations of reciprocal lattice 
vectors. ‘The extended zone scheme for the tin structure is extremely 
complicated and the very tedious remapping procedure was in part 
carried out for us by Petersen (private communication). Recently 
Harrison (1959) has described a much simpler procedure by which sections 
through the Fermi surface in the reduced zone scheme can be obtained 
directly and we have used Harrison’s method to construct the surfaces 
shown in fig. 5, from sections in the central (100), (110) and (001) planes. 
The surfaces have been drawn more or less to scale, except that sharp 
corners have been rounded off, both because such a rounding is to be 


expected on account of band gaps, and in order to make the diagram 
clearer. 


ee EE eee 

} The zone has been incorrectly described by Mott and Jones (1936); for a 
further discussion see Chambers (1956) and Barron and Fischer (1959). Space- 
group considerations of the tin structure have been given by Mase (1959). 
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The free-electron Fermi surface in the reduced-zone scheme. The principal 
dimensions of the zone are: 
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where a=5:82 & and c=3°18 A. 
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The predictions of the free-electron model are then as follows. The 
first zone is entirely filled and the second and third ones nearly so, except 
for two small pockets of holes centred on the points W in the second zone and 
a multiply-connected surface containing holes in the third zone. In the 
fourth zone the Fermi surface is one of two sheets and consists of a singly- 
connected surface 4(b) formed by the intersection of two ‘pancakes’ 
containing electrons; this surface is surrounded by a region of unoccupied 
states which is in turn bounded by the multiply-connected surface 4 (q). 
In fig. 6 we show the central (001) section through the fourth zone, 


The central (001) section through the fourth zone showing the relation of the 


surfaces in zones 4 (a) and 4 (b) (fig. 5) to one another. The shaded 
regions are occupied. 


from which the relationship of the two surfaces to one another should 
become clear. The surfaces in the fifth zone all contain electrons and 
consist of a long ellipsoid along the tetragonal axis, ‘pear-shaped’ 
surfaces centred on the zone corners H (two surfaces per zone) and 
surfaces similar to that in zone 4(b), but tilted and centred on the points 
Vv (four surfaces per zone). Finally we have in the sixth zone small 
ellipsoidal surfaces centred on the points V and still smaller ellipsoids 


centred on the points K, and inclined as shown; both sets of ellipsoids 
contain electrons. 
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We have indicated with broken curves in fig. 5 the principal orbits of 
maximum or minimum area of cross-section which would be expected to 
give rise to a de Haas—van Alphen effect. These orbits are labelled with 
small Greek letters and the associated de Haas—van Alphen periods from 
eqn. (1) are given in table 2. As Harrison (1960) has pointed out, his 
method of constructing the Fermi surface in the reduced-zone scheme also 
enables one to find the free-electron orbit masses (eqn. (2)) in a straight- 
forward manner, and these masses are also given in table 2. 


Table 2. The free electron model 


Orientation Zone, orbit iP(G=") 


[001] 6 x 10-7} 
4-5 x 10-7 
1:8 x 10-7 
2-6 x 10-8 
1-0 x 10-8+ 
0-87 x 10-8 


6:6 x 10-7 
5-5 x 10-7 
8-0 x 10-8 


Bo LO 


1-8 x 10-8 
0-75 x 10-8 


ve, de cd 


5-2 x 10-8 
1-4 10-8 


0-94 x 10-8 
28° from [100] : 
in (100) plane 6x 10-7¢ 


= 30° from [110] 
in (110) plane 3p 6 x 10-7+ 


+ Estimated values. 


§ 4. CoMPARISON WITH EXPERIMENT 


Before carrying out a detailed comparison with experiment we note 
that the number of periods to be expected on the basis of the free-electron 
model is about nine when the magnetic field is along the [001] direction, 
and considerably more at other orientations. In practice the maximum 
number found at any orientation is about seven (figs. 1, 2 and 3), suggesting 
either that the actual Fermi surface may not be quite as complex as the 
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literal free-electron model would predict or that some periods are being 
missed because of low amplitude. One can get some idea about which 
of the free-electron periods might be absent from considerations of the 
x-ray structure factor, which vanishes across the vertical {110} planes of 
the zone but not across the inclined {101} planes (Mott and Jones 1936, 
p. 164). This leads us to suggest that the (001) section through I of 
the actual Fermi surface might not be too different from that in the free- 
electron model, but small portions of the surface near the upper and lower 
ends of the zone might be considerably affected by the presence of energy 
gaps. In particular, the zones one to five are degenerate at the zone 
corner H in the free-electron model, but in going from H in the first zone 
to the corresponding point in the fifth zone we must cross four planes of 
the type {101} with non-vanishing structure factor; this means that we 
might expect the ‘pear-shaped’ surfaces centred on H to be reduced in 
size. Such considerations lead us to suggest that the small pockets of 
holes at W in the second zone may be filled, and that the two sets of 
ellipsoidal surfaces in the sixth zone together with the thin ‘intersecting 
pancake’ surfaces in the fifth zone may all be absent. Our simplifying 
proposals find some justification in the fact that we have found no 
experimental evidence for the particular period anisotropies expected 
for these surfaces; thus, for example, we should expect one of the larger 
ellipsoids in the sixth zone to exhibit a period maximum when the field 
is inclined at 28° to [100] in the (100) plane (i.e. field along the axis of an 
ellipsoid), but as can be seen from fig. 1 this does not happen. Again, 
the proposed simplifications result in a considerable reduction in the 
total surface area and so are, broadly speaking, in keeping with Chambers’ 
(1952) findings that the total surface area is only 0-43 times the free 
electron value. In what follows we shall therefore compare the results 
of figs. 1, 2 and 3 with the surfaces in only the four lower diagrams of fig. 5, 
but excluding the thin ‘intersecting pancake’ surfaces in the fifth zone. 

It is convenient to begin by considering first of all the long-period 
oscillations (upper halves of figs. 1, 2 and 3) which arise from the smaller 
portions of the Fermi surface. The period variations in both (100) and 
(110) planes of the groups A and B can be accounted for by surfaces of 
revolution parallel to the [001] axis; the angular variation of group A is 
approximately as the cosine of the angle measured from the [001] axis, 
thus suggesting a more or less cylindrical surface with its axis along 
[001], while that of group B is less rapid than the cosine suggesting a 
convex surface. According to the free-electron model we could therefore 
associate the A and B oscillations with the orbits § and pu respectively 
and we see from tables 1 and 2 that the observed periods are in fact 
reasonably close to the free-electron values. We do not expect to be 
able to follow the A and B oscillations all the way to the [100] and [1 10] 
orientations since the orbits 5 are then no longer defined while the orbit 
é' of the fifth-zone ellipsoid is not particularly favourable for a de Haas—van 
Alphen effect. These considerations suggest then that the C oscillations, 
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which are strongest at [100], are not connected with the groups A and B 
and we tentatively ascribe them to the orbits y of the third-zone surface. 
However the C oscillations extend all the way to [110] (C’ in fig. 3) and 
this is only possible if the slanting arms of the multiply-connected surface 
are not present (i.e. “pinched off’ near W). This proposal is supported 
by the fact that fig. 2 reveals no oscillations which can reasonably be 
ascribed to the B-orbits round these arms; the existence of B-orbits would 
require a period which had a maximum value at about 30° from [110] in 
a (110) plane (fig. 2). Finally, the very weak oscillations D are too vague 
to suggest anything definite. 

We now turn to the short-period oscillations (lower halves of figs. 1, 
2 and 3) which are related to much larger portions of the Fermi surface. 
The magnitudes of the three periods E, F and G observed at [001] agree 
reasonably well with the three short periods 7, c, ¢, respectively, predicted 
by the free-electron model. The G oscillations are fairly isotropic away 
from [001] and this might be accounted for by the fact that the free- 
electron orbit 7 in zone 4 (a) is actually concave where it crosses the orbit 
¢, so that the area of the orbit € would not be expected to vary much when 
the field is tilted some 30° to 40° from [001]. 

The short periods found near [100] and [110] have not, however, found 
a plausible interpretation in terms of the free-electron model. Thus the 
M oscillations, which probably have a maximum period at [100], might 
be associated with the orbits 7 in zone 4(a), but their amplitude variation 
(fig. 4) finds no plausible explanation. Again the H and H’ branches 
might be associated with the orbits € in zone 4(b), but they do not appear 
to join smoothly with the E branch as would be required by this inter- 
pretation. The results of fig. 1 and fig. 4 in fact strongly suggest that it 
is the L oscillations which join smoothly to E (ie. that it is L which 
must be associated with €) but the period minimum L at [100] does not 
fit in with the model. It can be seen that the free-electron model predicts 
three short periods near [100] corresponding to the orbits €, y and 7, 
while actually four, viz. H, K, L and M, are observed. Finally we 
should expect strong oscillations corresponding to the orbits «’ of zone 
4 (a) with a period maximum at [110], but no such maximum is observed ; 
the weak oscillations N, observed only at [110], have about the right 
period, but this may be fortuitous. The further weak period O also 
found only at [110] might be due to the orbit ¢ in zone 3. 

It is natural to enquire if slight modifications of the free-electron model 
might not help to eliminate the major discrepancies pointed out above. 
It turns out that only small changes are required to produce quite marked 
changes in the connectivity of the Fermi surface in the fourth zone. For 
instance the modification (i) of fig. 6 is rather attractive since it would 
remove the need for explaining the apparent lack of oscillations from the 
orbit ¢’, which would require a period maximum at [110]; although the 
orbit ¢ is then round a region of unoccupied states rather than occupied 
states, the area enclosed by the orbit remains practically unaltered by 
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this modification. On the other hand, modification (ii) would make the 
orbit 7 larger than the orbit ¢, requiring a still shorter period at [001] 
and two period maxima at [110], while (iii) would join orbits + and ¢ into 
one much smaller orbit which cannot now account for the two shortest 
periods at [001]. None of these modifications can, however, account for 
the presence of the four periods H, K, L and M, near [100] nor for the 
cylinder or hyperboloid of revolution about [100] which is required to 
explain the oscillations H, H’. In fact the only way in which such a 
cylinder or hyperboloid might fit into the framework of the free-electron 
scheme is by supposing the surfaces shown at the top of the fifth zone 
(fig. 5) to be reduced in size and yet to be so merged together as to form 
a multiply-connected surface of tubes whose axes lie along [100] axes 
through V. These [100] tubes could then give rise to the oscillations H 
while the oscillations F might conceivably be associated with a hole orbit 
around four tubes. This would also require an additional short period at 
[110]; the period L’ is of the right order of magnitude but it should be a 
maximum rather than a minimum at [110]. 

We conclude that the general degree of complexity of the experimental 
results is much the same as that of the free-electron model and that some 
of the details, notably the observed periods near [001], agree reasonably 
well with it. The serious discrepancies at other orientations, however, 
suggest considerable departures from the free-electron scheme, but it is 
hardly worth while speculating further about conceivable departures 
until more accurate and detailed data are available or until a better 
theoretical model is proposed with which the data can be compared. 
A ‘tight-binding’ calculation might be of help in suggesting surfaces to 
account for the short-period results in particular, and such a calculation 
is at present being carried out by Dr. M. Miasek. 


§5, Discussion OF OTHER RELEVANT EXPERIMENTS 


A certain amount of information about the Fermi surface in tin is 
available from other types of experiment and this will now be briefly 
reviewed. 

The anomalous skin effect was studied by Fawcett (1955) (see also 
Chambers 1956) and though it is now clear that the 3-ellipsoid model he 
used to interpret his results is far too simple, the interpretation of the 
present results is far from definite enough to allow detailed comparison 
with Faweett’s data. It is however worth mentioning that the free- 
electron surface of fig. 5 is quite capable of reproducing the salient 
features of the anisotropy in the surface conductivity. Thus the large 
radii of curvature of the orbits ¢ and y in zone 4 (a), especially in the region 
where they cross, would result in a high surface conductivity when either 
the [001] or one of the [100] axes lies in the surface of the sample, which 
is just what Fawcett finds. It is impossible at this stage to deduce a value 
of the total surface area for comparison with Chambers’ (1952) result of 
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0-43 times the free-electron value}, and the same holds for a comparison 
with the electronic specific heat. 

The most detailed study so far of cyclotron resonance in tin is that of 
Khaikin (1959) and, as can be seen from table 1, Khaikin’s masses at 
[001] are in fair agreement with the de Haas—van Alphen ones. However, 
the agreement is poorer at [100], but it is not clear if [100] is the diad 
axis to which Khaikin refers. The high value of m*=8-5m, found by 
Khaikin at an intermediate orientation is indicative of a large closed 
orbit extending through several cells of the reciprocal lattice and the 
free-electron model would suggest that this is most likely to be associated 
with the multiply-connected surface of zone 4(a). Assuming a free- 
electron-like interpretation of the present results, the de Haas—-van Alphen 
effective masses are in as reasonable agreement with the free-electron ones 
as one might expect (tables 1 and 2), bearing in mind the discrepancies 
observed in other metals (Harrison 1960). Cyclotron resonance experi- 
ments using circularly polarized radiation would be most helpful in 
deciding which of the surfaces contain electrons and which holes. 

The magneto-resistance results of Alekseevski and Gaidukov (1959) 
suggest that the Fermi surface of tin has open directions along [100] and 
[110] but this conclusion is based on only preliminary results. Further 
evidence on open orbits comes from ultrasonic attenuation experiments. 
Olsen and Morse (1959) find open orbits but no details are given. 
Mackintosh (private communication) also finds evidence for open direc- 
tions along [100] and [110]. The suggestion of open orbits along [100] 
is not inconsistent with the de Haas—van Alphen data (see pp. 1101, 1102) 
but we have found no period maxima at [100]; however, further comparison 
is premature until more detailed results are available. 
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Note added in proof.—We wish to point out that the free-electron surface 
of fig. 5 has been constructed assuming that electrons are Bragg-reflected at 
all proper Brillouin-zone planes whether or not the structure factor is zero ; 
a vanishing structure factor does not necessarily imply an energy gap 
which is identically zero over the whole plane (see refs. at foot of p. 1096 
Crystal symmetry requires the energy gap on the {110} zone faces, which 
have vanishing structure factor, to be zero only along the lines XL and XP 


+ Dr. Fawcett has informed us that his figure of 0-39 is low on account of a 
calibration error. 
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(Mase 1959; we are grateful to Dr. J. M. Keller for pointing this out to us). 
It appears that the only important alteration in the free-electron model 
brought about by this requirement is that an electron initially moving 
round orbit ¢ in zone 4 (a) will continue round half the orbit 6 in zone 3 after 
contacting the line XL of the (110) face, and so on. This results in a 
‘four-cornered rosette’ orbit normal to [001] whose maximum linear 
extension (along [110]) is somewhat larger than that of the zone but whose 
area is practically the same as the original orbit ¢. However such a 
‘freak’ orbit would have vanishingly small statistical weight and we have 
found no convincing experimental evidence for it. 


REFERENCES 


ALEKSEEVSKI, N. E., and GarpuKov, YU., P., 1959, J. exp. theor. Phys., U.S.S.R., 
37, 672: 

Barron, T. H. K., and Fiscuer, G., 1959, Phil. Mag., 4, 826. 

Brruincourt, T. G., 1952, Phys. Rev., 88, 242. 

CHAMBERS, R. G., 1952, Proc. roy. Soc. A, 215, 481; 1956, Canad. J. Phys., 34, 
1395. 

Conen, M. H., and Hutz, V., 1960 (to be published). 

Fawcett, E., 1955, Proc. roy. Soc. A, 232, 519. 

Gop, A. V., 1958, Phil. Trans. A, 251, 85. 

Harrison, W. A., 1959, Phys. Rev., 116, 555; 1960, [bid., 118, 1190. 

Heine, V., 1957, Proc. roy. Soc. A, 240, 340. 

Kuatrxin, M. S., 1959, J. exp. theor. Phys., U.S.S.R., 37, 1473. 

Masz, S., 1959, J. phys. Soc., Japan, 14, 1538. 

Mort, N. F., and Jones, H., 1936, Theory of Metals and Alloys (Oxford: 
Clarendon Press). 

OtseEn, T., and Morsg, R. W., 1959, Bull. Amer. phys. Soc. IL, 4, 167. 

SHOENBERG, D., 1952, Phil. Trans. A, 245, 1; 1953, Physica, 19, 791; 1957, 
Progress in Low Temperature Physics, Vol. IT, ed. C. J. Gorter (Amster- 
dam: North Holland Publ. Co.); 1960, Proc. Intl. Conf. on Fermi 
Surfaces in Metals (New York: Wiley). 


1105. =] 


The Determination of the Flux of Cosmic Ray 
Protons with Nuclear Emulsionst 


By C. J. WappiIneron 
H. H. Wills Physical Laboratory, University of Bristol 


[Received June 1, 1960] 


ABSTRAOT 


The flux of protons in the primary cosmic radiation has been determined 
in a stack of nuclear emulsions exposed during 1954 over northern Italy. 
This determination, which depended on the observation of the nuclear inter- 
actions produced by the protons, gave a flux at the top of the atmosphere of 
500+ 70 protons/m? ster sec, a value in good agreement with that expected 
from counter measurements. A simple model of the variation of the proton 
flux with depth in the atmosphere has been derived which appears to be in 
reasonable agreement with the experimental data. 


§ 1. INTRODUCTION 


Nuc.ear photographic emulsions have been extensively used to investigate 
the heavier components of the primary cosmic radiation, but have not 
provided much information on the primary protonst. Present data on 
these particles come entirely from measurements made with electronic 
counters. In particular, the Cerenkov-scintillator array of McDonald 
(1956) (see also McDonald and Webber (1959)) is almost entirely responsible 
for our present knowledge of the fluxes of these particles, since only this 
equipment appears able to discriminate adequately against the large 
background of upward moving albedo particles. Even this equipment, 
however, is unable to eliminate the re-entrant albedo, since these particles, 
which are predominantly electrons (McDonald and Webber 1959), have 
the same direction of flight as do those of primary origin. An alternative 
method of measuring these fluxes, and one which preferably discriminates 
against the re-entrant albedo, is thus very desirable. Such a method is 
provided by determining the flux from the nuclear interactions produced 
in nuclear emulsions. 
§ 2. MerHop 

When emulsions are used to investigate the heavier particles of the 
primary radiation it is usual to detect these nuclei by searching a given area 
for their tracks and from the frequency determining the flux. Such a 
procedure is hardly practical in the case of the primary protons due to 


+ Communicated by the Author. 
In what follows no distinction will be made between protons and the other 


isotopes of hydrogen that may be present in the primary cosmic radiation. 
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the extremely tenuous nature of the tracks left by fast singly-charged 
particles as they traverse the emulsion. Furthermore, even if by very 
careful scanning it were possible to detect all the tracks of particles 
producing minimum ionization which crossed a given area in the emulsions, 
it would be extremely difficult to distinguish the tracks left by protons 
from those left by other energetic but secondary singly-charged particles. 
Nor, in general, would it be possible to obtain any indication of the 
direction of motion of these particles. 

The nuclei of the primary radiation may, alternatively, be detected by 
observing the nuclear interactions, or ‘stars’, that they produce while 
passing through the emulsions. Until recently the determination of 
absolute fluxes from the frequency of these stars was also impracticable, 
due to the extreme difficulty of locating ali the interactions that are 
produced in any given volume. However, now that systematic observa- 
tions with beams of artificially accelerated protons are available, these 
detection inefficiencies can be calculated and allowed for. 

The additional requirement that the observed particles should produce 
a nuclear interaction has the great advantage that it serves to remove all 
those weakly interacting secondary particles, i.e. .-mesons and electrons, 
from the sample being considered, and thus greatly reduces the problem 
of how to take account of the albedo particles. The short lifetime of 
m-mesons means that the majority of those produced at any distance 
from the emulsions will have decayed to weakly interacting particles 
before they arrive. Thus only locally produced mesonic secondary 
particles are detected, and if the volume searched for stars is chosen to 
be near the top of the emulsion stack, then these will be predominantly 
upward moving, and can usually be eliminated from a study of the 
characteristics of the stars. The only serious contamination by secondary 
particles would appear to be that due to knock-on protons and to protons 
produced by the break-up of heavy nuclei above the emulsions. It appears 
to be possible to make a reasonable estimate of the abundance of these 
particles, and thus to calculate the flux at the top of the atmosphere. 


§ 3. EXPERIMENTAL DETAILS 


While no emulsion stacks which had been exposed at the same time as a 
Cerenkov-scintillor array were available in this laboratory, there was 
one which had a closely similar w-particle flux to that measured by counters. 
This stack was exposed on September 14, 1954 over northern Italy, and 
has previously been used to measure the a-particle flux and energy 
spectrum (Fowler and Waddington 1956), and the flux and charge 
spectrum of the heavier nuclei (Waddington 1957 a, Hillier and Rajopadhye 
1958). This stack consists of 80, 6 in. x 8 in. 600 4 G5 emulsions, and was 
exposed for 5-75 hours at an altitude of about 106000 ft. The flight 
curve has been published previously, e.g. the E stack (Rajopadhye and 
Waddington 1958). 
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A volume, all of which was less than 1 cm from the top of this stack, 
was scanned for stars which satisfied the following conditions. 


(i) The star lay within the inner 80° of the emulsion thickness, since 
Stars near an emulsion face are difficult to analyse. 

(ii) Each star had to have three or more black or grey tracks (N,,>3), 
one of which had a potential range more than 60 microns. This requirement 
on length removes those stars due to radioactive decays in the emulsion. 


Each of the stars found was carefully examined under high magnification 
by one or two observers for the presence of minimum tracks produced by 
fast primary or secondary particles. In those cases where such a track was 
found to lie within 60° of the vertical the star was, with a few exceptions 
discussed later, denoted as one produced by an energetic downward moving 
particle. The flux of these particles, J(X), which are assumed to be ail 
protons, that is, neglecting any downward meson component, is then 
given by 

J(X)= SOAS particles/m? ster sec a yore CL) 

V0)..4 

where A is the interaction mean free path of protons in emulsion, measured 
in metres, dN is the number of stars produced by these protons, V is the 
scanned volume in cubic metres, { is the solid angle over which particles 
are accepted and ¢ the duration of the exposure in seconds. This equation 
is derived assuming that the maximum distance a particle may traverse 
in the volume is much less than A. 

The value of dN obtained directly from scanning must be corrected for 
several forms of detection inefficiency. 


3.1. Detection Efficiency of Stars 


No attempt has been made to detect those stars with NV, <2, since it 
is well known that area-scanning cannot hope to detect such stars 
efficiently. Instead a correction has been made by using the experimental 
data available on the interactions produced in nuclear emulsions by 
artificially accelerated protons. Winzeler ef al. (1960) have published 
a N,, distribution for those stars found by ‘along the track scanning’ of 
6-2Bev protons. This distribution has been used to correct that found 
in this experiment. Figure 1 shows these distributions normalized for 
those stars with N,,>5, since it appears that there has also been some 
loss of stars with V,,=3 and 4 in this experiment. The correction factor 
derived for the loss of stars with NV, <5 is 1:56+0-13. 

The validity of this procedure depends on the assumption that the NV, 
distribution is nearly independent of energy, which appears to be generally 
true for energies above 1 or 2Bev. Furthermore, since the geomagnetic 
cut-off energy in this experiment was approximately 3-6 Bev for the 
protons, the median energy of these particles was about 5-7Bev. Thus 


4G2 
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the majority of the primary protons had energies which are not too 
different from those of the artificially accelerated protons with which 


they are being compared. 
Bigs i 


This experiment 


Normalized 6-2 Bev proton 


NO. OF STARS 


6 8 Le) 12 
NUMBER OF BLACK AND GREY PRONGS, N, 

The JN, distribution observed in this experiment and that observed by Winzeler 
et al. (1960) on 6-2 Bev protons. These distributions are normalized for 
those stars with V,>5. 

Fig. 2 
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The zenith angle distributions of those tracks which could definitely be assumed 
to be due to primary particles. The data has been divided into four 
equal frequency intervals. The mean of the first three, together with its 
standard deviations are shown by the horizontal lines. 
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3.2. Detection Efficiency for Primary Particles 


The tracks left by fast singly-charged particles are the most tenuous 
observable in nuclear emulsions and the most careful inspection is necessary 
to detect all such tracks associated with each star. Approximately 
one-third of all the stars were examined by two observers, and the good 
agreement between them suggests that both had very high efficiencies. 
However, such a comparison assumes that all tracks have an equal 
intrinsic probability of detection, and this may not be true. An alternative 
indication of the detection efficiency was supplied by studying the observed 
angular distribution of those particles which were assumed to be primary 
particles. It is reasonable to assume that the detection efficiency might 
be related to the zenith angle of the track, i.e. a particle dipping steeply 
in the emulsion might be more difficult to detect than one lying in the 
plane of the emulsion. The distribution of the space angles observed 
in this experiment is shown in fig. 2. It can be seen that if it is assumed 
that the incident particles are isotropically distributed there is a definite 
indication that particles with large space angles have been missed. To 
a first approximation these particles should be isotropically distributed 
and so a correction must be made for those missed. In this experiment, 
from a sample of 92 primary particles an additional 13-3 + 4-7 were assumed 
to have been missed. 


3.3. Identification of Primary Particles 


When a star has more than one minimum track associated with it, it 
is not always clear which was produced by the primary particle. Thus it 
is necessary to consider each star in detail in an attempt to determine 
which track has the greatest probability of being due to the primary. 

In an interaction of an energetic nucleon with an emulsion nucleus, 
the fast secondary particles are strongly collimated in the forward direction. 
As a result, in any particular star, the primary track may be taken as 
that which results in the least sum of the angles that the other minimum 
tracks make with the direction of motion of the assumed primary. In 
the majority of stars with three or more associated minimum tracks this 
condition is sufficient to determine the primary with a high degree of 
probability. However, in some cases no one track is clearly defined as 
being due to the primary. 

For this reason those stars having a minimum track within the solid 
angle of acceptance, 2, have been separated into seven groups. Sketch 
drawings of typical members of each group are shown in fig. 3. 


(a) Stars with just one minimum track. Since any albedo of energetic 
neutral particles can be neglected, this minimum track has a high 
probability of being that due to the primary. Ina sample of 100 stars, 
18% were of this type. 
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(b) Stars with two minimum tracks, one of which is in the lower hemi- 
sphere. While either track could have been produced by the incident 
particle, there appears to be a much greater probability that the one in 
the upper hemisphere should be that due to the primary (see later). 
Frequency 19%. 

(c) Stars with more than two minimum tracks, in which the only one 
in the upper hemisphere is clearly that due to the primary. Frequency 
Se 

(d) As (c), but with two or more tracks in the upper hemisphere, only 
one of which is clearly due to the primary. Frequency 8%. 

(e) As (d) but with no clear indication as to which track is due to the 
primary. These stars were given a weight depending on the number of 
tracks in the upper hemisphere within and outside Q. Frequency 8%. 

(f) Stars with more than two minimum tracks, in which one in the 
lower hemisphere is clearly produced by the primary. Frequency 5%. 

(g) Stars with two or more minimum tracks, all in the upper hemisphere. 
These were given a weight as for (e). Frequency 5%. 


Fig. 3 


Sketch drawings of typical stars. Each star has N,=4. 


Comparing groups (c) and (d) with group (f), it appears that only 
about 10% of the stars in groups (6) should have been produced by 
upward moving secondary particles. On a small sample of this sort 
this is a negligible fraction. It could, of course, be further reduced i 
exposing a smaller stack of emulsions, so as to reduce the amount of local 
matter available to generate these particles. 
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§ 4. EXPERIMENTAL RESULTS 


In scanning a volume of 1:47 cm?, 101 stars were found which had 
primaries in the solid angle of acceptance and V,>5. To these were 
added 14-6 + 5-1 stars to allow for the inefficiency of detecting primaries 
at large angles. The resulting 115-6 stars were multiplied by 1-56 + 0-13, 
the correction factor for those stars with V,<4. Then, using eqn. (1), 
with A=38-2+1-5em, the value quoted by Winzeler e¢ al. (1960), 
Q=7 steradians, and t= 2-07 x 104 sec, the time spent at ceiling altitude, 
J(X)=720 + 105 protons/m* ster sec, where J(X) was the flux of primary 
protons in the emulsion, and X = 10 g/cm? of residual matter. 


§ 5. GRowtH CURVE OF PROTONS IN THE ATMOSPHERE 


In order to extrapolate the flux value observed in the emulsions to the 
top of the atmosphere it is necessary to establish a growth curve for these 
particles. Furthermore, since in this experiment the emulsions were not 
re-orientated after they reached ceiling altitude it is essential to make a 
correction for those particles which entered while the emulsions were 
ascending and descending. For this purpose the growth curves must 
be known to appreciable depths in the atmosphere. It may be noticed 
that it is unnecessary to make a correction for those protons entering 
at sea level, since these will be closely collimated to the vertical, and 
it is unlikely that the stack would have been correctly orientated 
for them. 

There appears to be little or no experimental evidence regarding the 
flux of energetic protons as a function of atmospheric depth, and it is 
therefore necessary to construct a growth curve from a theoretical model. 
The simple, and somewhat artificial one chosen here assumes that in a 
typical interaction of an energetic proton with an air nucleus the proton 
is generally not absorbed but merely degraded in energy. ‘The energetic 
products of such an interaction consist mainly of 7-mesons, which on the 
average take about one-third of the incoming energy (Zhdanov et al. 1960) 
and occasional knock-on protons. Thus the proton multiplicity, n, is 
greater than unity for a number of successive interactions. Eventually 
the protons are either absorbed, or their energies are reduced below that 
necessary to produce tracks of minimum ionization. 

With these assumptions, if J,,(7, 0) is the flux of primary protons and 
protons produced by primary protons at a vertical depth x and zenith 
angle 0, and if J, is the flux of primary protons at the top of the atmosphere, 


then 
Hyo  Tiyhoh 


Tole, 0)=Jgexp (—ady)| 1+ + Be aa =JIoflx) (2) 


where A, is the interaction mean free path of protons in air, taken to be 
100 g/em?, and 7, 7, ... are the multiplicities in the first, second, ... 


interactions. 
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In this experiment, where particles were accepted at zenith angles of 
up to 60°, the flux at any depth « is given by 
iss J olY: 0) a 
vw & y 


[Zot 0) dw = 2x 


where y=2/cos 0 (Rossi 1952) : 
= 270 | “ehh dy =7J 


x 


pp 


which may be solved graphically. 


Fig. 4 
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The variation of the ratio Jpp/J with 2. 


A plot of the ratio J,,,/J9 against # is shown in fig. 4 for various cases. 
(a), (0) and (c) assume a three-stage process, i.e. use the first three terms 
of eqn. (2), with ny=n,=1-5; ny=n.=1-0 and n,=1-5, n»=1-0, respec- 
tively. (d) and (e) assume a two-stage process, with n,;=1-5 and n,=1-0 
respectively. 

It has been assumed that a three-stage process is more appropriate 
for protons above a cut-off energy of about 3-6 Bev, whereas a two-stage 
process applies for protons of lower energy. The only experimental 
justification for this appears to be the approximately correct prediction 
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of the sea-level proton flux mentioned later. A more refined model would 
have to take account of the range of energies of the primary protons. 
The values of n are somewhat more justifiable, since 1-5 and 1-0 represent 
reasonable limits and the combination of these two values in assumption (c) 
appears to represent the best model for this experiment. 

Approximately 50° of the nucleons in the primary cosmic radiation 
do not enter the atmosphere as free protons, but bound in nuclei. The 
median energy of these bound nucleons is approximately half that of the 
free ones, due to the influence of the geomagnetic field. For this reason 
it has been assumed that the protons released from these nuclei are capable 
of producing one fewer interaction than are the free protons before they 
are removed. For simplicity these bound nuclei have been assumed to 
enter only in the form of a-particles, and the primary «a-particle flux has 
been taken to be one-quarter of the primary proton flux. 

Thus, corresponding to a two-stage process for free protons 


J (x, 0) =2S x9 r i= [exp (—2/A,,) —exp (—2/A,)]=JS aofe(X) 


Pp 


where A, is the interaction mean free path of «-particles in air, assumed to 
be half A,, and J,, is the assumed flux of «-particles at the top of the 
atmosphere. 
So, as before 
2% 
Td ay = 20x | EEC) dy 
2 Y 
where J,,, is the flux of protons produced by a-particles in a two-stage 
process. 
Similarly, in a three-stage process 


A 
Jalt0) =2I 52 | exp (—x,) (<"- -1) 


Pp a 


So, as before 
2x 
Td yp = 202 i Fao) dy. 
z y 


Thus the total flux of protons, J(z), at any depth x is given by 
SQ) = J gat dap: 


Curves of the ratios J,,,/J49 and J(x)/J, as a function of x are shown in 
figs. 5 and 6. 

A partial test of the validity of this model is provided by comparing the 
flux predicted at sea level, i.e. 10A,, with that measured experimentally. 
Mylroi and Wilson (1951) have measured a flux of 0-55 protons/m? ster sec 
at Manchester, England. This may be compared with the value of about 
0-25 protons/m? ster sec predicted by models (d) and (e), which, because of 
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the low cut-off energy over Manchester, should be the appropriate models. 
In view of the simplified nature of these models this agreement appears 


to be remarkably good. 
Using the curves in fig. 6 it is possible to calculate the correction factor 


f(t) for the ascent and descent of the balloon, as well as that for the overlying 
atmosphere during the level flight, f(x). Values of f(¢) and f(x) are shown 
in table 1, together with the resulting values of Jo, given by 


Ty=JS(X). ft) - f(2); 


where J(X) is the flux measured in the emulsions (eqn. (1)). 


Fig. 5 
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The variation of the ratio J ap/ ao With a Curves a and b are for a 
three-stage process with n,;=1:5 and 1-0 respectively, while curve c 
is for a two-stage process. 


Thus, J(X) when corrected to the top of the atmosphere becomes, under 
the most plausible assumptions, model (c). 


J )=500 +70 protons /m? ster sec 
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where the quoted error does not include the uncertainties in the 
extrapolation to the top of the atmosphere. 
This value of Jy) may be compared with that of 


590 + 50 protons/m? ster sect 


Fig. 6 
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The variation of the ratio J(x)/J) with x. 


Table 1. Correction Factors 


ll 
f (x) F(t) -f (x) (protons/m? ster sec) 


485 +70 
0-753 540 +75 
0-694 500 + 70 
0-769 555 + 80 
580+ 85 


a EES Tar ErE 
+ After correction for the re-entrant albedo of about 40+10 particles/m? 
ster sec. 
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found by McDonald and Webber (1959) in a flight where the a-particles 
flux was almost precisely similar to that measured in this stack. 


§ 6. CONCLUSIONS 


The agreement between these two values of the proton flux appears 
reasonable, and suggests that there are no serious errors in either methods 
of determining these fluxes. 

The fluxes of all the nucleonic components of the primary cosmic 
radiation have now been measured in the stack of emulsions used in this 
experiment. These flux values are shown in table 2. 


Table 2. Fluxes of Nucleonic Components Measured in 
One Emulsion Stack 


Flux 


: Reference 
(particles/m? ster sec) 


Component 


oO 


This work 
Waddington (19 
Waddington (19 


Protons 
a-particles 
3<Z<5nuclei 
6<Z<9nuclei 
Z>10 nuclei 


i) 
5 


b) 
a) 


“ 
i 
7 
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With this demonstration of the feasibility of measuring proton fluxes 
with nuclear emulsions it should now be possible to use available 
experimental material to investigate the behaviour of this component 
of the primary cosmic radiation in some detail. However, it must be 
noted that emulsion experiments are still not able to give any information 
on the individual energies of these protons. Also that the accuracy of this 
method is limited by the uncertainties of the correction for secondary 
protons. This can be reduced by using small stacks exposed at great 
altitudes which have been re-orientated after reaching ceiling altitude. 
Under these conditions, and with the better data from machine experiments 
that should soon be available, accuracies of better than + 10% should be 
easily obtainable. 
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ABSTRACT 


A study was made of the intensity of scattering of monochromatic x-rays 
at small angles from a copper-base solid solution containing l5at. % 
aluminium. The intensity of scattering from a rolled foil was found to 
increase after heating at 150—250°c, a temperature range in which the alloy 
hardens, and there is little or no recrystallization. This increase is attributed 
to segregation of aluminium to the ribbons of stacking fault between partial 
dislocations. Various confirmatory tests were performed to check the 
validity of this conclusion. 


§ 1. INTRODUCTION 
1.1. Anneal Hardening 


THE present experiment was designed as an adjunct to a study (Chinowsky 
and Cahn 1961) of the anneal hardening of a-aluminium bronze. The 
cold-worked alloy becomes harder ifit is annealed below the recrystallization 
temperature ; a«-brass and some other solid solutions behave similarly. 
Hasiguti (1955) proposed the term ‘anneal hardening’ to denote this 
phenomenon. ‘Two alternative theories have been put forward to account 
for it : 

(1) The alloy may possess short-range order which is partly destroyed 
by deformation, and re-established by annealing. Jf the locking is in 
part determined by short-range order, this process can account for anneal 
hardening. 

(2) The dislocations formed during deformation may become locked 
by a segregation of solute atoms to the broad ribbons of stacking fault 
that separate the dissociated partial dislocations in an alloy of low stacking- 
fault energy. Such segregation has been postulated by Suzuki (1952). 
It can come about during an anneal following deformation, and cause 
the alloy to become harder. 

It is known (Borie 1958, private communication, Houska and Averbach 
1959, Davies and Cahn 1961) that «-Cu/Al] is liable to extensive short-range 
order and that this decreases with increase of the annealing temperature. 
Yet a sample quenched from 600°c, and accordingly almost disordered, 
does not harden on tempering at 200°C (a treatment which re-establishes 
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the order). For this, and other reasons, explanation (2) is more attractive, 
especially since «-brass and «-Cu/Al are known to have very low stacking- 
fault energies. More detailed discussion of the interpretation of anneal 
hardening will be found in a forthcoming paper by Chinowsky and Cahn. 


1.2. Small-angle X-ray Scattering due to Local Solute Segregation 


If ‘ Suzuki segregation’ does occur, it should be possible to detect it 
by means of small-angle scattering of x-rays. When a solid solution is 
freshly cold-worked, the new dislocations will be devoid of segregate. 
On subsequent annealing below the recrystallization temperature, solute 
will go to the stacking faults, which thus come to resemble GP zones in 
an age-hardened alloy. The local perturbations in the initially uniform 
electron density, resulting from a local concentration of solute, will in 
both cases give rise to an additional scattering of x-rays at small angles 
(Guinier and Fournet 1955). In the case of age-hardening alloys, such 
as Al/Ag and Al/Cu, this prediction has been amply confirmed (e.g. Walker 
and Guinier 1953). It was the purpose of the present experiment to look 
for this effect in an anneal-hardening alloy. 

The following conditions must be met if such an experiment is to have 
any hope of success : 


(1) The alloy must have solvent and solute of widely differing atomic 
numbers Z, since for a given dislocation density and a given degree of 
segregation, the intensity of small-angle scattering will be proportional 
to (fsolvent—Jfsolute)*, Where the f’s are atomic scattering factors. For 
small scattering angles, this = (Zsoivent — Zsolute)?- 


(2) The alloy must have a high dislocation density so that there are 
enough scattering centres for the predicted effect to be detectable. 


(3) The alloy must have a low specific stacking-fault energy, so that 
the ribbons of stacking fault between partial dislocations shall be wide; 
this again will enhance the scattered intensity. 


A copper-base alloy containing 15at.% aluminium, cold-rolled to 
50° reduction or more, satisfies all these conditions and was chosen as 
experimental material. The low stacking-fault energy and consequently 
large total area of stacking-fault ribbon in filings of Cu/Al alloys was 
confirmed by Paterson’s x-ray line shift method (Christian and Spread- 
borough 1957, Smallman and Westmacott 1957). Moreover, an analysis 
of the x-ray line broadening from Cu/Al filings, due to Guseva and Babareko 
(1959) yielded results closely similar to those obtained by Warren and 
Warekois (1955) in analysing line-widths from o-brass filings, and which 
they were able to interpret in terms of a faulted domain structure. 

Since the present work was completed, Nakajima (1959) has published 
a note in which he claims to have detected an increment of small-angle 
X-ray scattering due to Suzuki segregation in a Cu/Ni alloy. It will be 
surprising if this claim can be substantiated, since the alloy used by 
Nakajima signally fails to satisfy condition (1), above. 
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1.3. Double Bragg Reflection Mechanism 


A complication arises from the small-angle x-ray scattering which 
is observed directly after deformation, and which is predominantly due 
to repeated Bragg reflections from pairs of subgrains in the deformed 
structure (Beeman and Webb 1959, Neynaber et al. 1959, Fricke and 
Gerold 1958, 1959, Atkinson 1958). However, by careful comparative 
measurements of scattered intensity before and after heat treatment, 
this effect can be allowed for, and this has been done in the present experi- 
ments. 

§ 2. ExpeRIMENTAL Meruops 


The small-angle scattering apparatus used is shown diagrammatically 
infig.1. A beam of copper Ka radiation, monochromated by an asymmet- 
rically cut bent quartz crystal, passed through the alloy foil before being 


Fig. 1 
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Schematic diagram of small-angle scattering apparatus. 


blocked by a beam stop at the focus. The alloy foil was set in a holder so 
arranged that the same area could be irradiated after the foil had been 
temporarily removed for heat treatment. The intensity of the beam 
was continuously monitored by a Geiger counter which measured the 
intensity of the (311) line diffracted by a very thin cold-rolled foil of 
aluminium, set in the beam behind the quartz crystal. 

The intensity scattered by alloy foil was measured by means of a pro- 
portional counter (20th Century Electronics Ltd.) set behind a slit on 
a rotating arm. Used in conjunction with a pulse amplifier and a pulse 
height analyser, this counter had a cosmic-ray background of only 4 
counts/min. This rather elaborate counting equipment was essential 
since the scattered intensity was so small that a typical Geiger counter 
background of about 1 count/sec would have been of the same order as 
the ‘signal’. Fluctuations of the beam intensity were allowed for by 
accumulating a fixed number of monitor counts, and counting the pulses 
arriving at the proportional counter during the corresponding period. 


P.M. 4H 
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The high purity alloy used, containing 15 at. % aluminium, was et 
and carefully homogenized, reduced by further cold rolling to a ie 
100 thickness, annealed and rolled down to 50y, and finally thinne 
uniformly to a thickness of about 15 by means of chemical polishingy. 


§ 3. RESULTS 


Figure 2 shows a scattering curve obtained from a cold-rolled alloy 
foil, using in front of the proportional counter a narrow slit that subtended 
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Small-angle scattering curve from an alloy foil rolled to 50 °%% reduction. Correc- 
tion for background count and parasitic scattering has been made. 


10’ of arc at the sample. The results in fig. 2 have been corrected for 
background count and for residual parasitic scattering, which was 
appreciable for scattering angles less than 1:8°. This curve is typical 


} The polishing bath, used at 70°c, contained 50% glacial acetic acid, 30% 


oe nitric acid, 10° cone. hydrochloric acid and 10% cone. orthophosphoric 
acid, 
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of those obtained for a rather weakly deformed metal (Neynaber ef al. 1959, 
Webb and Beeman 1959). 

The sample was then annealed progressively at 250°o and the scattered 
intensities at several angles between 1° and 1-7° were measured after 
each anneal. The intensities progressively increased, until after 13 hours’ 
annealing all had increased by about 20°. This increase was statistically 
significant in spite of the fairly large probable error in counting. 

Further experiments were now made using a larger slit at the counter, 
subtending 30’ of arc at the sample, and keeping to a fixed scattering 
angle of 2°, measured at the slit centre. In this way a higher statistical 
accuracy could be achieved. The intensities quoted hereafter are uncor- 
rected for parasitic scattering (which was very weak at 2° scattering angle). 

A foil was reduced 70% in thickness by rolling, cut in half and chemically 
thinned, and the intensity of x-ray scattering from each piece was measured. 
In each case this was about 470 counts during the period needed by the 
monitor to accumulate 30000 counts. Most of this intensity was presum- 
ably due to double Bragg scattering. The pieces were then annealed in 
stages, one at 200°C and the other at 250°c, and the scattering from the 
same areas of the foils measured after each anneal. Figure 3 shows the 
results. The points denote the increase in the scattered intensity in 
terms of counts per 30000 monitor counts ; the vertical line through 
each point represents the probable statistical error. The solid curves 
represent results for the isothermal hardening of a similar massive alloy 
sampley, plotted with a suitably normalized ordinate scale so that the 
hardness and intensity plots should concide as nearly as possible. 

Figure 3 shows the following : (1) The intensity scattered from both 
samples increases significantly, the maximum increase from both samples 
being about 15%. (2) There is a very rapid initial rise in the scattering 
power ; an anneal of less than one minute at 200°c suffices to produce 
an 8% increase in intensity. (3) At 250°C, recrystallization begins after 
some hours for the degree of deformation used, and the scattered intensity 
then falls steeply. (4) The hardness increase at each stage is proportional 
to the increase of scattering power. 

Attempts were made, using the point-focusing small-angle camera 
described by Atkinson (1958) (kindly loaned by Dr. P. B. Hirsch) to record 
the small-angle scattering photographically so as to study its distribution 
in azimuth, but these attempts failed because of the low intensity available. 


§ 4. Discussion 
The results shown in fig. 3 might in principle be associated with any 
of the following changes in the specimen during anneal hardening : 
(1) The re-establishment of short-range order. 
(2) Growth of subgrains in the alloy and subsequent modification 
of the intensity of double Bragg reflection. 
(3) Suzuki segregation of solute to the stacking faults. 
+ Results due to Mr. 8. Chinowsky. 
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Alternative (1) is improbable since increase of order will diminish the 
weak Laue monotonic scattering (Warren and Averbach 1953) and should 
thus actually lead to a slight decrease of small angle scattering. As a 


Fig. 3 
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Isothermal kinetics, at 200 and 250°c, of the change in scattered x-ray intensity 
at 2°, and of the Vickers pyramid hardness. 

X-ray results: Cu/15 at. % Al alloy, rolled 70%. Ordinate shows increase in 
counts per 30 000 monitor counts, above initial value of 470 counts. 
Hardness results: Cu/14 at. % Al alloy, rolled 75%. Ordinate shows increase 
in Vickers hardness number. 


check, a piece of the alloy foil was annealed at 600°c and quenched, which 
recrystallized and disordered it, and the scattered intensity at 2° was 
measured. The foil was then annealed 17 hours at 150°c, which suffices 
to produce substantial short-range order (Davies and Cahn 1961). The 
scattering intensity after this treatment had not changed, within the 
statistical counting error of +3%. Alternative (1) is thus ruled out. 
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Alternative (2) can also be excluded fairly definitively. The theory 
of double Bragg reflection (Webb and Beeman 1959) shows that the total 
scattered flux due to the double reflection mechanism is proportional to 
the inverse square of the mean subgrain diameter ; subgrain growth should 
thus lead to a rapid reduction of intensity scattered at small angles. This 
has indeed been confirmed by Reizman (1959), who studied the change, 
on annealing, of small-angle scattering from rolled copper foil. The 
resultant subgrain growth led to a steady decrease of the small-angle 
scattering ; there was no sign of a transient increase. We carried out 
an experiment to confirm this. Copper and alloy foils, rolled to approxi- 
mately the same reduction, were each annealed 17 hours at 150°c, which 
is enough to harden the alloy considerably but not to cause any recrystal- 
lization in either the alloy or the copper. The x-ray scattering at 2° 
from the unannealed alloy was 927 + 20 counts per 63 000 monitor counts, 
and it was 1076+ 20 counts after heat treatment—an increase of 16%. 
The corresponding counts for the copper were 791 and 784 respectively ; 
the difference is not significant. A further anneal of 20 hours at 250°c 
caused substantial recrystallization of the copper and the count fell to 
625. There is no reason why the double Bragg reflection process should 
be differently affected by heat treatment in the case of the alloy than in 
the case of copper. 

There is therefore no cause, theoretical or experimental, to suppose 
that subgrain growth might in itself cause an increase of the small-angle 
scattering from the alloy. 


§ 5. CONCLUSION 


Having eliminated the alternatives, we are left with the conclusion that 
the increased small-angle scattering from the cold worked alloy, resulting 
from heat treatment at 150—250°c, is due to a segregation of solute to 
the stacking faults. 

The present results are merely indicative; more detailed exploration 
is desirable of the isothermal kinetics of the change of the scattered intensity, 
of its dependence on the degree of deformation, on the state of polarization 
of the scattered x-rays, etc., but the low x-ray intensity available from our 
equipment precluded such a systematic study. In particular, experiments 
with rolled single crystals could show whether the distribution in azimuth 
of the scattered intensity was consistent with the presence of lamellar 
segregates parallel to the octahedral slip planes. 
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ABSTRACT 


The precipitation kinetics of a gold 1:5%, by weight, cobalt alloy have 
been followed by measuring the saturation magnetization of the precipitated 
cobalt and the results are compared with current theories. Marked differences 
were found between the ageing curves of air and water-quenched specimens. 
The remanent magnetization, coercive force and field to reduce remanence 
to zero of aged specimens have been measured and the results are discussed 
in terms of the single magnetic domain model. The magnetic properties are 
compared with those of the gold nickel system and the conclusion is reached 
that the dominant contribution to the anisotropy in the gold cobalt system is 
due to particle shape. In an appendix the remanence curve and coercivity 
for a random array of single domain ferromagnetic particles having a linear 
distribution of anisotropies have been calculated and shown to predict the main 
features of the observed results. 


§ 1. INTRODUCTION 


THE magnetic properties of fine ferromagnetic precipitates in non- 
ferromagnetic matrices have been studied extensively in recent years and 
interpreted in terms of single domain theory. Much attention has been 
paid to the copper-cobalt system (Bate et al. 1955, Becker 1958, Mitui 
1958), where the cobalt-rich precipitate is thought to form as the face- 
centred cubic allotrope (Rodbell 1958). Unfortunately the analogous 
system, copper-—nickel, shows complete miscibility of the two constituents 
so that no comparison between the behaviour of the two face-centred 
cubic ferromagnetic metals in a copper matrix is possible. However, the 
gold-nickel and gold-cobalt systems both show limited solid solubility. 
The lattice parameters of the face-centred cubic forms of cobalt and nickel 
are very similar (3°545kX and 3-517kX), thus leading to similar dis- 
registries between precipitate and gold matrix (lattice parameter 4-704 kX), 
of about 12%. The spontaneous magnetizations of cobalt and nickel do, 
however, differ considerably (164 and 58 e.m.u. per g), and any difference 
in magnetic properties between the two systems will be in large measure 
due to this. In particular, the theory of Stoner and Wohlfarth (1948) 
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predicts that if the coercivity of an array of small single domain particles 
is due to shape anisotropy then it will be proportional to the spontaneous 
magnetization of the material. On the other hand, contributions due to 
crystalline and magnetostriction anisotropy will be inversely proportional 
to the spontaneous magnetization. A comparison of the two systems, 
with different spontaneous magnetizations but similar metallographic 
properties, should therefore indicate the dominant contribution to the 
observed anisotropy. The gold-nickel system has been investigated by 
Lothian e¢ al. (1958), and results on the gold—cobalt system are presented 
in the present paper. 

According to the phase diagram of Raub and Walter (1950) it is possible 
to retain 1:5°%, by weight, of cobalt in solid solution in gold at 960°c, and 
then to form a ferromagnetic precipitate, containing about 99% of cobalt, 
by ageing at a lower temperature. The allotropic transition from the 
face-centred cubic to the hexagonal close-packed form takes place at a 
temperature less than 500°c in cobalt-gold alloys containing 99% cobalt 
(Koster and Horn 1952). If, therefore, ageing is carried out at 500°C or 
above, the equilibrium precipitate will have a face-centred cubic structure. 


§ 2. PREPARATION OF SPECIMENS 


A gold—cobalt alloy containing 1-52%, by weight, of cobalt was made by 
melting the constituents together in an argon arc furnace on a water-cooled 
copper hearth. The alloy was then cast into a rod, about 0-4cem in 
diameter and 5 cm long, by melting in vacuo in an evacuated silica tube 
using an induction furnace. This rod was homogenized, in vacuo, for 
36 hours at 960°c before disc-shaped specimens about 0-1 em thick were 
cut from the rod. 

The specimens were annealed for 5 min at 960°c in evacuated silica 
capsules, and quenched to room temperature before ageing at 565°c. 
Some specimens were water-quenched, by withdrawing the silica capsules 
rapidly from the furnace and breaking them under cold water; others 


were quenched more slowly by allowing them to cool in air without breaking 
the tube. 


§ 3. MEASUREMENTS OF SATURATION MAGNETIZATION 


Magnetization at fields greater than 40000e was measured using a 
Sucksmith ring balance (Sucksmith 1929 and 1939). 

The magnetization of a quenched unaged specimen was proportional 
to field at room temperature but showed a tendency towards saturation 
at liquid nitrogen temperature. Quenched specimens showed no remanent 
magnetization or coercive force. 

After ageing for 3min the magnetization curve did not extrapolate 
through the origin, there was a remanent magnetization, and a coercive 
force of 9000e was required to reduce the magnetization to zero. 
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In order to determine the saturation moment of the ferromagnetic 
precipitate a correction was made for the paramagnetism of the matrix 
solid solution by subtracting the small magnetization of the quenched 
specimen from that of the aged ferromagnetic specimen at each measuring 
field. The saturation magnetization was then determined from a linear 
extrapolation against the reciprocal of the measuring field. 

The increase in saturation magnetization with ageing time for water- and 
air-quenched alloys is shown in fig. 1. The maximum value of the 
saturation moment was 1:72+0-02e.m.u./g which compares with the 
1-8+0-le.m.u./g predicted for this alloy from the data of Raub and 
Walter (1950). 

Fig. 1 
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Variation of saturation magnetization with ageing time. 


The initial rate of precipitation in water-quenched specimens was much 
greater than in air-quenched specimens. This is almost certainly due to 
there being a greater number of nucleation sites in the water-quenched 
alloys where the faster quench will retain the vacancies produced at the 


annealing temperature of 960°C. 


§ 4. Remanenr MAGNETIZATION 


The remanent magnetization after the application of a field of 17 000 0e 
was measured using a torsion balance as described by Bate et al. (1955). 
The ratio of the remanent magnetization to the saturation magnetization 
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was determined for 29 specimens including both air- and water-quenched 
specimens. There was no significant change in this ratio with ageing 
times from 1:5 to 50 min, the mean value of the ratio being 0-52, the 
standard deviation among the specimens being 0-05. The ratio of 
remanent to saturation magnetization is expected to be § for a random 
array of single domain precipitate particles having predominantly uniaxial 
anisotropy (Wolhfarth 1955). 

Complete remanence curves were plotted for a few specimens using the 
torsion balance. After ageing the specimen possesses zero magnetization, 
the individual particle moments lying along easy directions of magnetiza- 
tion. On applying and removing a small field, H, the magnetization of 
particles of low anisotropy is reversed and the specimen is left with a 
remanent magnetization J,(H). Further increase of H increases I,(H) 
until a saturation remanence, J,(00), is reached where all the particles 
have a component of magnetization parallel to the direction of the applied 
field. Ifa reverse field is now applied the remanent magnetization, J ,(#), 
decreases to zero at H,. and to saturation remanence in the reverse direction 
at still higher fields. 


Fig. 2 
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A complete remanence curve for a water-quenched specimen, aged for 
7 min, is shown in fig. 2. A field of more than 10 000 0e was required to 
produce saturation remanence indicating that the Specimen must have 
contained particles with intrinsic coercivities as high as 1000008. 
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(‘Intrinsic coercivity’ of a particle is used to indicate the coercivity the 
particle would exhibit if the direction of lowest magnetic anisotropy 


Solid line - theoretical prediction 
X — experimental observations 
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Fig. 4 
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Variation of coercivity and remanence coercivity with ageing time for 
air-quenched specimens. 


energy were parallel to the field.) It was of interest to compare these 
results with an expression of Wohlfarth (1958) which gives: 


I4(H)=1,( 0) —21,(H). 
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This formula is easily proved for any array of non-interacting single 
domain particles of uniaxial anisotropy. In fig. 3, J4(#H) is plotted against 
I,(H) showing a marked deviation from the ideal. The deviation may be 
due to the fact that the particles are not strictly uniaxial but possess 
mixed uniaxial and cubic anisotropies. 


Fig. 5 
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Variation of coercivity and remanence coercivity with ageing time for 
water-quenched specimens. 


§ 5. ComrctvE Force AND FIELD TO REDUCE REMANENCE TO ZERO 


Coercive forces greater than 900 oe were measured using the Sucksmith 
ring balance; forces less than 9000e were determined with the torsion 
balance. The field to reduce remanence to zero, H,, was also determined 
with the aid of the torsion balance. The variation of the two fields with 
ageing time in air-quenched specimens is shown in fig. 4. The maximum 
coercive force of 2050+ 500e was attained after about 10 min ageing; 
the maximum H,, of 28000e being reached after about the same time. 
The ratio of the two fields was thus about 1-36 after 10 min ageing. For 
a random array of uniaxial single domain particles of uniform anisotropy 
this ratio should be 1-09 (Wohlfarth 1958); a ratio of 1-36 indicating a 
range of particle anisotropy of about the same order as observed by 
Lothian et al. (1958) in Au-Ni alloys. The change of the fields with time 
is consistent with the growth of single domain particles whose shape 
anisotropy increases as their size increases. The fall in the fields after 
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the optimum ageing time is reached is then due to the onset of multi- 
domain behaviour when the particles reach a critical size. Water-quenched 
specimens behave similarly, as shown in fig. 5, reaching a maximum H P 
of 2050 oe after about 10 min ageing but attaining a higher value of H , at 
35000e. The ratio H,/H, is thus about 1-7, indicating a greater range of 
particle anisotropy than is found in air-quenched specimens. 

The coercivity of the aged specimens was increased by c. 6000e on 
cooling to liquid nitrogen temperature. A graph of H, against temperature 
for a water-quenched specimen aged for 10 min is shown in fig. 6. 
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Variation of coercivity with temperature for a water-quenched specimen 
aged for 10 min. 


§ 6. ErrEctT oF STRAIN 


A preliminary experiment indicated that the aged alloys behaved in a 
similar way to Au—Ni alloys (Lothian et al. 1958); that is, the coercivity 
is greatest perpendicular to the axis in a drawing experiment and parallel 
to the axis in a compression experiment. Thus the coercivity of an 
air-quenched Au—Co specimen increased from 1880 oe to 2200 oe parallel 
to the axis and decreased to 15100e perpendicular to the axis after a 
compressive strain of 6%. This contrasts with the observations of 
Livingston and Becker (1958) in a Cu—Co alloy where they observed 
enhanced coercivity parallel to a tension axis, and the increase in anisotropy 
energy produced by the strain was consistent, up to 20% strain, with 
homogeneous deformation of matrix and precipitate. ‘The contrast in 
behaviour may be due to the high disregistry of 12% between precipitate 
and matrix in the Au-Ni and Au—Co systems compared with the 1-7% 
disregistry in the Cu-Co system: if the disregistry is small the particles 
will remain coherent with the matrix and dislocations will pass easily 
through the particles thus making shape deformation of the particles 
possible. If disregistry is large the particles quickly lose their coherence 
and the passage of dislocations is impeded by the particles. 
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§ 7, X-RAY MEASUREMENTS 


Debye-Scherrer photographs of quenched and aged needle specimens 
were taken using copper radiation. No extra lines due to the precipitate 
phase were detected. The quenched specimens gave sharp reflections at 
all Bragg angles, and the lattice parameter of the solid solution was 
4:0535kX. Ageing an air-quenched specimen for 2-75 min caused the 
high-angle reflections to blur toward the low Bragg angle side of the 
photograph; this blurring becoming more pronounced with further 
ageing. After 18min, precipitation was almost complete, and after 
2500 min ageing the high angle lines were again sharp, the matrix having 
a lattice parameter of 4:0658kX. The water-quenched specimens behaved 
similarly. The films show that as precipitate particles begin to form the 
matrix far from the particles is left unchanged in composition and a 
cobalt-depleted region surrounds the particle. This suggests that almost 
as soon as precipitation begins the full lattice disregistry of 12% must be 
accommodated at the matrix particle interface. As ageing proceeds the 
concentration gradient in the matrix becomes smaller until the cobalt 
concentration reaches its equilibrium value and the high angle diffraction 
lines become sharp. 


§ 8. Discusston 
8.1. Saturation Magnetization 


The saturation magnetization of the specimen is a direct measure of 
the amount of precipitate and it was of interest to compare the kinetics 


of the precipitation process with a formula of Zener (1949) which gives the 
fraction precipitated as : 


or[-(9" 


where + and m are constants and ¢ is the ageing time. Thus if 
log log (o../¢.—0,) is plotted against log t, where o,, is the saturation 
magnetization after complete ageing and og; is the saturation magnetization 
after time ¢, it should give a straight line of slope m. In fig. 7 this 
function is plotted for both air- and water-quenched specimens. The 
graphs show that the increase of saturation magnetization with time (in 
minutes), can be well represented up to 60% precipitation, in air-quenched 
specimens by: 

o,= 1-72[1—exp — (¢/7-50)?* 44] 
and in water-quenched specimens by 

o,= 1-72[1 —exp — (¢/1-66)*23}], 


Ham (1958) has treated the problem of precipitation from a supersaturated 
solid solution in a rigorous way and has shown that a formula of the type 
derived by Zener should be a good approximation for the first 50% of the 
precipitation process. He has further shown that the exponent m will 
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be 1:5 for any shape of particle in a diffusion-limited process, where surface 
and strain-energy effects are negligible, and in which the inter-particle 
separation is about ten times the particle diameter. The eccentricity of 
the particles in such a process will remain constant after a short initial 
transient. For situations where a surface reaction is dominant m tends 
to 3, while for diffusion-limited growth where the initial particle volume 
is greater than one-tenth the final volume, m will be between 1-0 and 1:5. 


Nig e7 


+2 


Logio logy) (G0 /O% — ot) 


oO os 1:0 1:5 2:0 


Log,, (time, min.) 


Logarithmic plot of saturation magnetization against ageing time. 


For the air-quenched specimens m>1-5, indicating that precipitation 
cannot be purely diffusion-limited and surface and strain-energy effects 
are probably important. In the water-quenched specimens m<1-5 so 
that the nature of the nuclei formed by rapid quenching must be such that 
the coherency strains are more easily accommodated and large particles 
are formed in the initial transient process. 
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8.2. Coercivity and Field to reduce Remanence to Zero 


The two fields H, and H, reached their maximum values after about 
10 min ageing in both air- and water-quenched specimens. The optimum 
value of H, was, however, significantly higher in water-quenched specimens 
although H, was approximately the same for air- and water-quenched 
specimens. The larger value of H,/H, in water-quenched specimens 
indicates that they contain a wider range of particle anisotropies than the 
air-quenched specimens. This must be associated with the fact that as 
the optimum ageing time is reached in water-quenched specimens the 
supersaturation is small, precipitation being about 93% complete, but in 
air-quenched specimens it is larger since precipitation is only 55% complete. 
According to Nabarro (1940) a low supersaturation favours increasing 
eccentricity of the growing particle. Thus the shape anisotropy and 
hence the intrinsic coercivity of particles growing from a solution of low 
supersaturation will reach higher limiting values before the size of the 
particles allows domain walls to form. The range of anisotropies will 
therefore be wider in water-quenched specimens since particles of high 
anisotropy are favoured. 

A graph of H, against temperature, for a water-quenched alloy aged 
for 10 min is shown in fig. 6. The variation of H, is linear over the range 
77°K-380°K. Stacey (1959 a) has derived a formula relating the magnetic 
energy of a single domain particle to its volume v, and the coercivity of a 


random array H,: 
H=vH I, +kT InCr 


where J, is the spontaneous magnetization per cm® of the particle, & is 
Boltzmann’s constant, 7' is the absolute temperature, 7 is the time taken 
to measure H,, and C is a frequency factor whose maximum value at room 
temperature is given as 6 x 10sec by Stacey (1959b). Thus assuming 
I, is independent of 7’ over the range considered : 


dH dH. 
(The factor k/7' arises because Coc 7’, Stacey 1949b.) If now the magnetic 
energy is due to shape anisotropy alone then dH#/dT is zero and the variation 
in coercivity arises only from temperature dependent fluctuations of the 
particle magnetization. Hence 


k(n Cr + 1/7) 
I,dH aT ~ 


From the graph dH,/d7' = — 2-87 oe/c deg and taking 7 as 1 sec (this is 
not critical), gives v=1-010 x 10-8cm? and hence a particle radius of 
order 624. This calculation gives an average size for the particles 
contributing to the observed coercivity. If the direction of low shape 
anisotropy is associated with a direction of low crystalline anisotropy 


Magnetic Study of Precipitation in a Gold—Cobalt Alloy 1137 


energy, the contribution to dH,/dT' due to cubic crystalline anisotropy 
will be —0-66 (Bean e¢ al. 1959). 
Thus 
= = —0-66vl,; 
therefore 
k(n Cr + 1/T) 
I,(dH,/dT + 0-66) ’ 


leading to a slightly larger value for v as 1-311 x 10-8cm3 and radius 
68 A. 

The air-quenched specimens show approximately the same slope so that 
the particle radii are not appreciably different in air- and water-quenched 
specimens aged for 10min. The value obtained above for the particle 
radius agrees with that found by Becker (1957) for the onset of ferromagnetic 
properties in copper—cobalt alloys. 


8.3. Remanent Magnetization 


The remanent magnetization of annealed specimens after the application 
of a field of 17 000 oe was half the saturation magnetization. This suggests 
that the precipitate particles have predominantly uniaxial magnetic 
anisotropy (Wohlfarth 1955). In fig. 2 the remanent magnetization was 
plotted against the applied field for a water-quenched specimen whose 
mean coercivity was 15700e and field to reduce remanence to zero was 
27000e. The saturation remanence was not reached until a field greater 
than 100000e was applied, showing that particles having intrinsic coer- 
civities of this order were present. In the appendix a remanence curve 
for a simple linear distribution of anisotropies is calculated and compared 
with the experimental curve and using the linear distribution giving 
the best fit with the experimental remanence curve a value of H, was 
calculated as 1440 oe. 

The saturation remanence for an air-quenched specimen, for which 
H,, was 2500 oe and H, 1850 oe, was reached at about 7000 oe, thus showing 
that a narrower range of anisotropies was present in these specimens. 

The intrinsic coercive force of a prolate spheroid is given by (V,—N,)Lo 
where NV, is the demagnetizing factor along the short axis and NV, along 
the long axis (Stoner and Wohlfarth 1948). The maximum value, for 
an infinitely long spheroid is thus given by 27J,. For cobalt this value 
is 8962 0e but a spheroid of axial ratio 10 reaches 98% of this maximum 
value. A maximum contribution to the particle coercivity of 2K/I 
is possible due to crystalline anisotropy. Ifthe precipitate is face-centred 
cubic this adds a maximum contribution of c.14000e making a total of 
c.104000e. The observed maximum intrinsic coercivity can therefore 
be explained in terms of a cubic precipitate showing no magnetostriction 
effects. A magnetostriction contribution to the coercivity is possible 
if the precipitate experiences a uniaxial stress. This will only occur 
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if the precipitate remains partially coherent with the matrix. In the gold— 
cobalt system the disregistry is about 12% so that the precipitate will 
begin to break away from the matrix when the particle diameter reaches 
100/12 atom spacings since a dislocation can then be introduced to relieve 
the strain between the two lattices. Thus the precipitate will begin 
to break away from the matrix when its radius becomes greater than 
304, which is about half the lower limit of particle size to show ferro- 
magnetic, rather than super-paramagnetic, behaviour (Becker 1957). 
The precipitate will therefore have started to break away after the optimum 
annealing time and there is not likely to be a large magnetostriction 
contribution to the magnetic anisotropy. 

As shown in fig. 6, the remanence curve shows a marked deviation from 
the relation proposed by Wohlfarth (1958), relating the remanence J,(H) 
gained by the initial application of a field to an aged specimen and the 
remanence J,(H) remaining after the application of a demagnetizing 
field to a previously saturated specimen. 

The relation will only be strictly obeyed for an array of perfectly uniaxial 
particles rather than the particles of mixed cubic and uniaxial anisotropy 
which are formed in practice. It is possible for particles in their initial 
state to have magnetization vectors not lying in the uniaxial direction 
but in one of the cubic easy directions. 

The application and removal of a suitable field will turn these vectors 
into the easier uniaxial directions. On applying a demagnetizing field 
the cubic easy directions may no longer be available as remanence 
positions since there is no suitable field which will bring the magnetization 
vector close enough to the cubic axis for it to rest there on removal of the 
field ; the critical angle for reversal of the vector in the uniaxial direction 
being greater than the critical angle for return to the cubic direction. 
This means that the demagnetizing mode of the remanence curve will 
be harder than the initial magnetizing mode so that the remanence coer- 
civity predicted from the initial magnetizing curve will be less than that 
found by demagnetization. 


8.4. Comparison with Au—Ni 


The maximum coercivity found by Lothian ef al. (1958) in aged Au—Ni 
alloys was 3250e. It was, however, observed that the coercivity increased 
if the alloy cross section was reduced by cold drawing. Thus Robinson 
(1957) found that after as little as 13% reduction the coercivity both parallel 
and perpendicular to the drawing axis increased to about 6000e. The 
coercivity perpendicular to the axis was, however, greater than that 
parallel to the axis and further reduction increases this difference, the 
coercivity parallel to the axis decreasing. The initial rapid increase in 
coercivity can be explained if it is supposed that the elongated particles 
are in tension, thus giving a negative contribution to the coercivity. 
A small amount of cold work is then sufficient to relieve this tension as 
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dislocations cross the particles. Orientation effects will, however, tend 
to make the coercivity parallel to the drawing axis less than that per- 
pendicular to the axis. 

The maximum coercivity observed in the Au—Co system was 2050 oe and 
in the compression experiment the coercivity parallel to the compression 
axis increased but perpendicular to the axis it decreased. An initial overall 
increase in coercivity, both parallel and perpendicular to the axis, was 
not therefore observed but the orientation effect was strong. Thus it 
appears that the magnetostriction contribution is much less than the 
shape anisotropy contribution in Au-—Co. It is to be expected that the 
magnetostriction contribution will be smaller for cobalt than for nickel 
alloys since it is inversely proportional to the spontaneous magnetization. 

The fact that the maximum observed coercivity in the cobalt alloys 
is bigger than in the nickel alloys suggests that the shape anisotropy term is 
dominant since this is proportional to the spontaneous magnetization. If 
it is further supposed that the particle shape is determined mainly by 
metallographic factors the constant of proportionality will be roughly 
the same for both systems because of the similarity of cobalt and nickel 
in structure and atomic size. With these assumptions it is possible to 
compare the two systems and estimate the magnetostriction contribution, 
to the anisotropy using the relation : 


= Al + = +2 
where J, is the spontaneous magnetization and A, B, and C are constants 
determining the shape, magnetostriction and magnetocrystalline contri- 
butions. Thus assuming A and B to be the same for both the nickel and 
cobalt systems and allowing 700o0e and 100 oe for the magnetocrystalline 
contributions in cobalt and nickel, gives : 


$ s 


1350 = 14004 + for cobalt, 


1400 


225= 500A+ a for nickel. 


This gives A = 1-04 and B= — 1-47 x 10°. 

The negative value of B is expected for nickel and cobalt in tension since 
the magnetostriction coefficient is negative for nickel in all directions and 
for cobalt in (111) directions. Using the formula of Stoner and Wohlfarth 
(1948) for the magnetostriction contribution from a random array of 


particles : 
3Ab 


$s 


H,=0-479 


where \ is the magnetostriction coefficient and 6 is the uniaxial tension, 
gives a value of about 3-0 x 10! dynes/cm? for 6. 

This is small compared with the upper limit for b suggested by Stoner 
and Wohlfarth of about 2 x 10” dynes/cm?. 


412 
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The value of A gives a mean value for the axial ratio of the elongated 
particles of approximately 1-6. 

A suggested breakdown of the anisotropy contributions for the two 
systems is shown in the table : 


H Shape | Crystalline Magnetostriction 


c 


Co 2050 1455 700 — 105 


Ni 325 520 100 —295 


Thus relaxation of tension, by slight cold working, in the Au—Ni alloys 
increases the coercivity to 6200e which agrees with experiment. In 
the Au—Co alloy the increase in coercivity, due to relaxation of tension 
by cold work, is swamped by the changes due to the orientation of the much 
larger shape and crystalline anisotropy contributions. 
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APPENDIX 


REMANENT MAGNETIZATION AND COERCIVITY FOR A LINEAR 
DISTRIBUTION OF PARTICLE ANISOTROPY 


Current theories based on the single magnetic domain hypothesis 
generally deal with arrays of particles of uniform anisotropy. In practice, 
and especially in heterogeneous alloys, a wide range of anisotropies is 
present and theories based on a single anisotropy must prove inadequate. 
The simplest distribution of anisotropies which can be considered, and 
whose parameters can be easily varied, is a linear distribution with an 
upper and lower cut-off value. Such an approximation is sufficiently 
good to indicate the main features of an array containing a range of 
anisotropies and in particular to show how the ratio H,/H, depends on this 
range. 

For a linear distribution of particle anisotropy the volume of precipitate, 
dv, whose intrinsic coercive force lies between Hy and (H,+dH,) is given 
by : 

dv=AdH, for values of Hy between H, and H, 
where A is a constant and H, and H, are the lower and upper cut-off 
values of Hy. Such a distribution, with suitably chosen values of H, 
and H,, gives quite a good approximation to the results obtained in the 


present work and shows how the ratio H,/H, varies with the range of aniso- 
tropies present. 
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Remanent Magnetization 


The remanence dJ,(H) for a random array of particles of intrinsic 
coercivity Hy can be tabulated from the data of Stoner and Wohlfarth 
(1948) using the method of Wobhlfarth (1955). This gives the remanence 
on first applying a field H to a completely random array and the field H, 
is that which produces a remanence of half the maximum value. 


Thus : 


dl,(H)= Af(H/H,)dH 
therefore 
al, (H 1\ 
( = Aftia(¢ ) 
where h=H/H, ; therefore 
1/h=H2/H 
a =4 | HOAL ae 
ve 1/h=Hy/H 
Fig. 8 
Ir (H) 
Ir (00) 
0-8 
0:6 
0-4 
O-2 
° ate 
0:5 1-0 5 2:0 h 


Variation of reduced initial remanence with reciprocal of reduced field for a 
random array of uniform anisotropy. 


In fig. 8, f(2) is plotted against 1/h and can be used to evaluate the integral 
of eqn. (1) for any particular value of H. When H=dH, all the 
particles have turned toward the field direction and eqn. (1) becomes: 


vi 1h=1 
eee | f(h)d(1[h) = A [1— H,/H). 
2 1/h=Hi/H 
Equation (1) can arin sa re-written as : 


a“ me" fi) ach) (2) 
[,( 0) aia aif (Rc na se 
The field H, is that for which I,(H)/I,(0©)=0-5 ; therefore 


0-5= 


= a {1-855 — (H,/H,)} for H,/H, <1, (H/H,)>2 
ee 1 


Hence H,=0-269 (H,+H,) for H,/H,<1, H,/H,> 2. 
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The values of H, obtained experimentally from the initial and demagnet- 
izing branches of the remanence curve were not the same ; the value 
from the demagnetizing branch being the greater. In fig. 9, therefore, 

Fig. 9 


Ia (CH) 
Ir (<0) 


0-8 


Solid curve - theory 
© - observed values 


-15 -10 
Field, kilo-oersted 


-0:4 


-O°6 


Comparison of demagnetizing remanence curve, for a linear anisotropy array, 
with observed values for a water-quenched specimen aged for 7 min. 


the experimental value of J4(H)/J,( 0c) from fig. 2 is compared with that 
for a linear array of anisotropies with H,=500 and H,=9530 using the 
relation of Wohlfarth (1958) : 

I,(H)=I,( co) — 21,(H) 
to calculate I;(H). The theoretical curve has been chosen to fit at H, 
and as closely as possible elsewhere. The uniform anisotropy distribution 
fits least well at the higher fields indicating that the true distribution 


contains fewer particles of anisotropy less than H, than is implied by the 
linear distribution. 


Coercivity for a Linear Array 


Stoner and Wohlfarth (1948) have tabulated the reduced magnetization 
curve for a random array of particles of intrinsic coercivity H 9. Thus 
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by a similar argument to the above, for a linear anistropy distribution : 
I(H) H 1h=H9/H 

= | G(hyd( lit) Pee & 10 (3) 
1/h=Ai/H 


I(co) H 2— fi, 
where J(#) is the intensity of magnetization for the demagnetizing branch 
of the curve. The reduced magnetization, g(), for a random array of 
particles of uniform anisotropy is plotted against 1/h. the reciprocal of 
the reduced field, in fig. 10. The coercive field is a function of H, and 


Fig. 10 


0-6 


0-2 


> pa 
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Variation of reduced magnetization, along demagnetizing branch, with 
reciprocal of reduced field for a random array of uniform anisotropy. 


H, and is that field for which the integrand of eqn. (3) is zero. It can be 
represented to an accuracy of better than 1% by : 


H,=0-133H,f{ 1+ 2-5(H,/H,)} ; 


further 
H,,=0-269 H,(1+H,/H,) fon ial, boli, > 2.4 
therefore 
14+ 4H,/H 
=p 1 2 . 
seh ek od ai + 2-5H,/H, 
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This formula is valid for values of H,/H, up to 0-37 when H,/H, becomes 
equaltoone. For values of H,/H, greater than one H,, has to be calculated 
using fig. 8 to determine the appropriate integrand. In fig. 11 the ratio 
H,/H, is plotted against all possible values of H,/H,. 


Fig. 11 


1-0 , He 
fo) Oo2 04 O06 O8 ro 0H 


Variation of H,/H, with ratio of minimum to maximum anisotropy field in a 
linear anisotropy distribution. 


For the water-quenched specimen considered above the observed values 
of H, and H, were 2700 and 15700e, so that H,/H, was 1:72. For the 
best fitted remanence curve the suggested value for H,/H, of 0-0524 
gives a value for H,/H, of 1-88 and for H, of 14400e. Better agreement 
could doubtless be obtained by adopting a more realistic anisotropy 
distribution but the present approximation is sufficiently good to show 
the main features to be expected from an array of ferromagnetic particles 
having a range of anisotropies. 

The calculations suggest that 2-02 will be the maximum possible value 
of H,/H, for such an array. Values of H,/H, greater than this could be 
produced in alloys containing a large proportion of super-paramagnetic 
material of zero coercivity. Thus almost all the values of H,/H, observed 
by Lothian et al. (1958), in Cu-Co alloys, were greater than 10 indicating 
a large super-paramagnetic contribution. However, in Au-Ni alloys 
and in the present work H,/H, was always less than 2, indicating that 
almost all the precipitate was in ferromagnetic form. 
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ABSTRACT 


The dynamic stress-strain behaviour of aluminium crystals fatigued in a 
push-pull machine at 13 cycles/second is outlined. Room temperature 
interruptions to tests carried out at —196°c, —78°o and 20°o resulted in a 
transitory hardening which is attributed to the movement of point defects 
to dislocations. 


A study has been made of the behaviour of 99-99% aluminium single 
crystals fatigued in a push—pull machine at 13 cycles/second. The stress 
pulse was fed to the X-plates of an oscilloscope, while the strain was mea- 
sured by a capacitance strain gauge connected to the Y-plates so that a 
dynamic hysteresis loop was obtained continuously during a fatigue test 
(fig. 1, Pl. 146). In this way the dynamic hardening of the crystals could 
be studied over a wide range of fatigue conditions. 

In tests carried out at room temperature, the width of the hysteresis 
loop (dynamic plastic strain) obtained at a given alternating stress gradually 
decreased (fig. 1 (a), (b) and (c)) in the way observed for much lower fre- 
quencies by Thompson ef al. (1955). The maximum hardening occurred 
within the first 5°% of the fatigue life, and a steady state was reached after 
about 15-20% (fig. 1(c)). However, in tests carried out at —78°c and 
—196°o, no cyclic hardening whatsoever was detected, the loop width 
remaining at a steady value for large numbers of cycles (> 200 000 cycles). 

In these low temperature tests, it was found that intermediate resting at 
room temperature resulted in a hardening of the crystals when fatigue was 
resumed at —196°c or —78°c. Figure 2 (a) (Pl. 146) shows the hysteresis 
loop of a crystal 22B after 6500 cycles (+ 1015 g/m? resolved shear stress) 
at —78°c; the crystal was then rested at room temperature for 60min 
and the loop width on immediate re-testing is shown in fig. 2(6). After 
a further 100 cycles (fig. 2 (c)) the loop width had reached a size comparable 
to that prior to the interruption. The hardening tended to be most 
pronounced in the initial stages of a test when the deformation was carried 
out at —196°c or —78°c (fig. 3), but a large effect was still obtained on 
resting a week after 4 x 104 cycles. 


+ Communicated by the Author. 
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In tests undertaken at room temperature, not only was the observed 
hardening greater (fig. 4(a)) but the effect persisted longer aE can be seen 
by comparing the behaviour of an identical crystal at —196°c (fig. 4 (0)). 


Fig. 3 
104 iB 
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1-2 +1015 gms/mm* RESOLVED SHEAR STRESS 
2 thr, 
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nase 10 20 30 40. x103 


CYCLES 


The effect of interruptions at room temperature on the loop width measured 
at —78°c. 


Fig. 4 
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° -196°C 
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(a) The effect of room temperature interruptions on loop width measured at 
room temperature, (b) the effect of room temperature interruptions on 
the loop width measured at —196°c. (a) and (6) relate to tests carried 
out on identically oriented crystals. 
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The purity of the metal used is sufficiently high to exclude the possibility 
of impurity atom locking as a cause of the observed hardening, and to 
indicate instead point defects. Broom and Ham (1957, 1959) have suggested 
that point defects contribute to fatigue hardening by locking dislocations 
and moreover have shown yield phenomena when pure copper crystals 
previously fatigued at —196°c were deformed in tension in the range 
—196°c to 20°c. This view is supported by the electrical resistivity 
measurements of Rosenberg (1958) on copper fatigued at —253°c; the 
observed increase in resistivity is large with respect to that obtained as a 
result of uni-directional deformation and anneals out from —100°c 
upwards. The present experiments also indicate that point defects are 
formed during fatigue and migrate during resting at room temperature. It 
is thought that vacancy movement to dislocations will also explain the 
fatigue hardening observed in the present tests at room temperature. 

The hardening obtained on resting is of a transitory nature for on further 
cycling the loop width (plastic strain) rapidly approached that which 
existed prior to resting. This appears to be in direct contrast to the 
results of tensile tests on fatigued crystals where a permanent upward 
displacement of the stress-strain curve was observed. The much bigger 
effect on the loop width observed at room temperature suggests that higher 
concentrations of point defects are generated than at lower temperatures, 
which would be in accord with the much heavier development of slip 
striae at room temperature (Hull 1958). 
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ABSTRACT 


The total flux and energy spectrum of primary cosmic-ray alpha particles 
at the top of the atmosphere have been measured in emulsion stacks flown 
during September 1956 and June 1958. Values of 240+26 and 130+15 
particles/m? steradian second respectively were obtained for the total flux. 
This reduction in the flux during June 1958 was caused by the partial 
removal of those alpha particles with energies between 100 and 2000 
Mey/nucleon. These changes have been correlated with the 1l-year cycle of 
solar activity. 


§ 1. INTRODUCTION 


It was pointed out by Forbush (1954) and others (e.g. Singer 1958), on the 
basis of observations made with neutron monitors at ground level, that 
periods of high solar activity, as measured by the Zurich sunspot number, 
were associated with a low cosmic-ray intensity. The corresponding 
variations in the flux and energy distribution of the primary radiation 
have been observed with nuclear plates flown to high altitudes and Freier 
et al. (1958) were able to show that the alpha particle flux in 1957 (measured 
over Minnesota) was reduced to about 50% of its value in 1954. Their 
results also indicated that this reduction was not caused by the application 
of an abrupt cut-off to the energy spectrum but that particles with energies 
as low as 200Mev/nucleon were still able to reach the emulsions. The 
object of this present work has been to obtain further data on the energy 
spectrum and total flux of primary alpha particles using balloon-borne 
nuclear emulsions. The results have been obtained during September 
1956 (Experiment 1) when the level of sunspot activity had risen consider- 
ably above the minimum of 1954 and in June 1958 (Experiment 2) near the 
maximum of the current solar cycle. 


§ 2. EXPERIMENTAL PROCEDURE 
2.1. September 1956. Hapervment 1 


An emulsion stack consisting of 150 plates, each 40 cm x 30cm x 600 p in 
volume was flown to a mean height of 117000ft over Minnesota (North 
America) at a geomagnetic latitude of 55°N. for 64 hours. At this height 
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the thickness of overlying air is equivalent to 6-5g/cm?. The observation 
and measurement of the alpha particles was divided into two parts, the 
first being the determination of the differential energy spectrum below 
250 Mev/nucleon and the second the measurement of the total flux and the 
energy spectrum up to 1000Mey/nucleon. For each section of the work a 
different sample of plates was searched using appropriate scanning criteria. 
In the first section only those tracks of density greater than 6 x minimum 
ionization were selected along a line 5mm from the top edge of each plate. 
In the second section this density criterion was extended to 3 x minimum 
in order to include alpha particles of all possible energies. The alpha par- 
ticles themselves were selected from all tracks recorded in the scanning by 
measuring the multiple coulomb scattering and grain density of each track. 
This method of identification was similar to that used by Fowler et al. 
(1958), except that in the present case the ionization was measured visually 
instead of using a densitometer. The scanning efficiency was checked in 
various ways. Some of the scanning was repeated by different observers 
and there was found to be no evidence for the omission of alpha particles. 
In addition, the number of alpha particles was found to be independent of 
the zenith angle (after correcting for atmospheric absorption) and the 
depth at which the track was found. It was assumed, therefore, that the 
scanning loss was negligible. 


Fig. 1 
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Blob-gap density—range relationship for alpha particles. 


The methods used to obtain the alpha particle energies differed in 
parts 1 and 2 of experiment 1. In part 1 the measurement was based 
on a determination of either the grain density or the range of each track. 
Kach alpha particle was followed through the emulsion stack until it either 
interacted, left the stack or came to rest. In the case of those particles 
which stopped, a measurement of the distance between the point of entry 
and the end point was sufficient to determine their energy. In those 
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cases where the track did not end, the blob-gap density (Fowler and 
Perkins 1955) was measured either in the last few millimetres of available 
track or near to the point of entry. The blob-gap density was then 
related to the particle energy in the following way. A number of 
calibration plates were selected from those already scanned in part 1. 
These plates were searched for relativistic alpha particles which were 
used to give a value for the minimum blob-gap density in the plate 
concerned. All those alpha particles of known residual range which 
traversed the calibration plates were then blob-gap counted in those 
plates, to obtain a relationship between the blob-gap density, expressed as a 
multiple of the minimum density, and the residual range. The result of 
this procedure is shown in fig. 1. The full line shown in the figure was used 
to give the residual range of those alpha particles for which a direct range 
measurement could not be made. The energies which were obtained by 
this method were probably accurate to about + 10%. 

In part 2, the alpha particle energies were obtained by a multiple 
coulomb scattering measurement, which was made over a length of track 
at least greater than 8mm, using 500, cell sizes. The momentum and 
hence the energy of the particle was calculated using a value of the 
scattering constant obtained from the calculations of Fichtel and 
Friedlander (1958) and by eliminating the random noise between 500 
and 1000p» cell sizes, it being assumed that the contribution of the noise 
to the observations varied linearly with the cell size. This method gave 
a standard deviation of about + 15% in the energy region considered. 


2.2. June 1958. Hxpervment 2 


The dimensions of the stack in this second experiment were 
20cm x 12cm x 600 and it was flown with the 20cm edge uppermost 
over the same place as in Experiment 1 for 8} hours. The height was such 
that the mean thickness of air above the emulsions was 4g/cm*. The 
plates were scanned for all tracks greater than 3x minimum, and the 
energies of the alpha particles were determined to an accuracy of 10-15% 
by multiple scattering measurements. 


§ 3. RESULTS 
The total fluxes were calculated using eqn. (1). 
—— 11ii1 
J=Nsecd.7 Oo 
J = Flux of alpha particles ; 
N =Number of particles ; 
A = Area of emulsion ; 
T =Time at maximum altitude ; 
Q = Solid angle ; 
sec 0 = Mean value of secant of zenith angle. 


particles/m? ster. sec. WAS Jatt Mia 9) 
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The values of V which were used in eqn. (1) were obtained from the numbers 
actually observed by correcting for the absorption of alpha particles and 
their production from the interactions of heavier nuclei in the emulsion 
and overlying atmosphere. These corrections were applied by inserting 
the fluxes of nuclei heavier than helium into the diffusion equations of 
Kaplon et al. (1954) and using the mean free paths and fragmentation 
probabilities of Rajopadhye and Waddington (1958). For the data of 
Experiment 1, values of 5, 15 and 5 particles/m? ster. sec were taken for the 
light, medium and heavy nuclei respectively, (cf. Dainton et al. (1952) 
and Kaplon et al. 1954), but for the results of Experiment 2 these values were 
reduced by 50%, it being assumed that the heavier components of the 
cosmic radiation decreased in the same proportion as the alpha particle 
flux (cf. Freier et al. 1959a, b). The correction for those particles which 
arrived during the ascent of the emulsions was eliminated in Experiment 1 
by launching the stack upside down and in Experiment 2 by launching it in 
an overturned position. The stacks were preset to rotate into the correct 
orientation when the maximum altitude was reached. Unfortunately, 
in Experiment 1, a short time was spent at altitude before rotation took 
place and a small correction was needed to allow for this. 


Experiment 
number 


Number of alpha Flux at top of 
particles observed atmosphere 


Date 


September 18, 1956 : 240 + 26 particles/m? 


ster. sec 
June 14, 1958 130+ 15 particles/m? 
ster. sec 


The flux values which were obtained are shown in the table in which the 
errors shown are the standard deviations on the number of tracks observed, 
together with an allowance for the uncertainty in the thickness of the 
emulsion at the time of exposure. The number tabulated for Experiment 1 
does not include the 75 particles with energies less than 250 Mev/nucleon 
which were obtained in part 1. 

The differential energy spectra were obtained by extrapolating the 
observed energy of each alpha particle to the top of the atmosphere and then 
dividing the particles amongst suitable energy groups. The fluxes for 
each energy interval were found using eqn. (1) with the appropriate 
corrections. The resulting energy spectra are shown in fig. 2. There is 
good agreement between the two sets of low-energy results obtained in 


Experiment 1 even though the energies of the particles were found by 
different methods. 


§ 4. Discusston 


The results for the total flux of alpha particles obtained in the present 
experiment are compared in fig. 3 with other values obtained at a similar 
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geomagnetic latitude but at different periods in the solar cycle (see 
Waddington (1960) for a list of references). The same figure also shows the 
variation with time of the Zurich sunspot number (Chernosky and Hagan 
1958). Itis clear that there is a strong anti-correlation between the alpha 
particle flux and the solar activity, as indicated by the sunspot number, so 
that periods of high solar activity are associated with a low flux of alpha 
particles. The most rapid change in the alpha particle flux coincides with 
the sharp rise to high sunspot number between 1955 and 1957. 


Fig. 2 
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Differential energy spectrum. 


The energy dependence of the variations can be seen from fig. 2. The 
spectrum obtained in September 1956 showed a pronounced peak in the 
region of 350Mev/nucleon. At energies higher than that value the experi- 
mental points appeared to follow the power law given by eqn. (2) and shown 
by the hyphenated line in the figure : 


dJ — Constant , 


See rae 2 
dE (1+£)?5 ’ (2) 
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E=alpha particle energy in units of the particle rest energy. te 
energies below the peak there was a rapid decrease in intensity as the 
particle energy decreased to 100Mev/nucleon. The alpha particle energy 
spectrum has been measured over a similar energy range In J une 1954 ata 
geomagnetic latitude of 60°N. and, over the same place as this experiment, 
in October 1950 (Fowler et al. (1957) and Fowler et al. 1958), and in August 


Fig. 3 
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1956 (McDonald 1958). The results of McDonald were obtained with a 
Cerenkov and a scintillation counter and are in good agreement both with 
the results of Fowler et al. in 1950 and with the present experiment. How- 
ever, the spectrum obtained in 1954 differs from those just mentioned in 
that it rises to a higher value at the peak energy of 350 Mev/nucleon and 


Total Flux and Energy Spectrum of Cosmic-ray Alpha Particles 1157 


decreases less steeply between 350 and 150 Mev/nucleon. This may be 
related to the very low sunspot activity of 1954 compared with the 
increased activity in 1950 and 1956 on either side of the minimum. The 
spectrum of 1956 also showed that there was an appreciable number of 
alpha particles reaching the earth with energies approaching the limit of 
about 100Mev/nucleon set by the air and packing materials above the 
emulsions. The fact that particles with energies as low as this were 
observed provided an upper limit to the cut-off energy over Minneapolis in 
1956 which agrees with that observed by Fowler e¢ al. (1957) in 1950 and 
which is somewhat lower than that given by the recent calculations of 
Quenby and Webber (1959). 

The peak at 350 Mev/nucleon was still present in the 1958 results but it was 
much less prominent and the energy spectrum began to depart from the 
power law of eqn. (2) at a much higher energy, as can be seen in fig. 2. 
These results agree with those obtained in 1957 (Freier e¢ al. 1959b and 
Waddington 1960) and show that even during the years of maximum 
solar activity there are still an appreciable number of primary particles 
reaching the earth with energies as low as 200 Mev/nucleon as was observed 
by Freier et al. (1958). The height at which the emulsions were flown was 
sufficient to detect alpha particles of energy as low as about 80 Mev/nucleon 
but the number of particles observed with an energy of less than 
200 Mev/nucleon was insufficient to determine the cut-off energy. 

A number of attempts have been made to explain these variations but 
none are entirely satisfactory. Parker (1958b) has attempted to explain 
the phenomena by invoking the solar wind hypothesis of Biermann 
(1951, 1952 and 1957). According to this theory, magnetic lines of force, 
emanating from the sun, are carried towards the outer regions of the solar 
system by a continual flow of ionized gas. It was predicted (Parker 1958 a) 
that this flow would become unstable in the region between the orbits of 
Mars and Jupiter and produce a shell of disordered magnetic field which 
would present a barrier to the incoming cosmic ray particles. Parker was 
able to show that for a given particle energy, the flux (J) in the vicinity of 
the earth would be related to that outside the barrier (J..) by eqn. 3. 


Jo(H)=JoBexp {—4.(5) -cae ash ree 


In this equation Z/M is the change to mass ratio of the particles con- 
cerned and A is a constant depending only on the magnetic and geomagnetic 
properties of the barrier. The results of calculations based on eqn. (3), and 
assuming that J,,(#) was given by eqn. (2), are shown in fig. 4. The 
lower curve was obtained using a value of A(Z/M)* of 0-58 which is that 
required by the calculations of Parker. ‘The upper curve uses a value for 
A(Z/M?) equal to 0-12. There is fair agreement between the theoretical 
predictions and the experimental results. Unfortunately the theory does 
not predict any variation in the energy spectrum above about 2 Bev/nucleon 
in contrast to the observations made using neutron monitors (Forbush 


1954, 1958). 
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Comparison between theory and experiment. 


§ 5. CONCLUSIONS 


The total flux and the energy spectrum below about 2000 Mev/nucleon 
of primary cosmic ray alpha particles have been measured on two separate 
occasions. ‘The results exhibit a strong anti-correlation between alpha 
particle intensity and solar activity and are in general agreement with 
the theory of Parker (1958b). In spite of this, the theoretical understand - 


ing of the relationship between the cosmic-ray observations and the solar 
activity is still unsatisfactory. 
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ABSTRACT 


The emission in the 40—-10004 range, excited in a potassium target by 
electron bombardment, has been examined using the grazing incidence photon 
counting spectrometer described previously. A total of seven lines and bands 
are observed of which two are identified as the L and M valence spectra, 
three are identified as diagram lines, and two broad lines remain unidentified.. 
An energy level scheme for solid potassium in emission has been constructed. 

A note regarding the energy of the L, level in sodium, magnesium and 
aluminium is appended. 


§ 1. INTRODUCTION 


THE M spectrum of solid potassium has been observed previously by 
Kingston (1951) and by Fisher ef al. (1958) while Tyrén (1937) has 
observed the Ll, Ly doublet. Although the present research has been 
concerned principally with valence band spectra, in the case of potassium 
a number of atomic lines also fall in the relevant spectral range and their 
wavelengths may be measured. These lines can be identified by reference 
to the extrapolated energy values for the L and M levels, and hence exact 
values for these energy levels may be obtained. 


§ 2. INSTRUMENTAL 


The spectrometer has been described in detail previously by Fisher eé al. 
(1958). A number of modifications have been made subsequently and are 
discussed elsewhere by Crisp and Williams (1960). They include a cold 
baffle which nearly surrounds the target and which practically eliminated 
carbonaceous contamination (Crisp 1958, Ennos 1953, 1954), an improved 
photomultiplier detector, and replacement of the aluminium grating by a 
ruled glass grating which gave intense spectra down to 45 A and enhanced 
reflectivity in high orders (Crisp 1960). 


+ Communicated by the Author. 
{ Now at National Research Council, Ottawa. 


1162 R. 8. Crisp on 


The targets were prepared by the evaporation of pure potassium, from 
a stainless-steel cup, on to a water-cooled copper target. The metal was 
loaded into the cup, freshly cut under the packing oil and wet with oil 
which was subsequently pumped off. In the preparation of an actual 
target, before depositing the specimen metal, the specimen chamber was 
‘oettered’ by a massive evaporation of potassium and the copper target 
freshly cleaned by means of the scraper. Pressures of lower than 
10-°mm Hg were produced in this way and cleanliness of the specimen 
was ensured. 

Wavelengths were calculated from the fiducial marker ‘pips’ on the 
chart records. The method of calculating the theoretical calibration 
curve has been indicated by Fisher et al. (1958) and Crisp (1960). 
A correction curve to this calibration was drawn using the emission 
edges of the L spectra of aluminium and magnesium, and the peak of 
the boron K spectrum in several orders. Skinner’s (1940) wavelengths 
were taken as the ‘true’ wavelengths of these features and unknown 
wavelengths could then be determined absolutely with a probable error 
of +0-5A4. Differences in wavelength may be measured with an error 
of +0-054. The random errors incurred in making actual measurements 
on the charts may generally be reduced to a negligible figure by averaging 
a large number of records. The small systematic error introduced by the 
finite width of the analyser slit is automatically eliminated in the 
construction of the correction curve. 

The slits were set at 40-45, and on the criterion of Fisher (1954) this 
corresponds to an optical resolution of 0-94. The instrumental width 
referred to hereafter is merely this 0-9 converted to electron volts and 
will be somewhat smaller than the true ‘line width’ as used for example 
by Skinner (1940). 


§ 3. EXPERIMENTAL RESULTS 


The seven observed lines could be conveniently divided into three 
groups: 


(i) The valence M spectrum (fig. 17). 


(ii) A number of diagram lines and the valence L spectrum (fig. 3). 


(iii) Two unidentified broad lines designated X and Y for convenience 
(fig. 3). 


The numerical data are set out in table 1, and the lines observed are 
illustrated in figs. 1 and 3 which are photographs of typical chart records. 
On each record the spectrum and order, target current and voltage, full 
scale counting rate and percentage error (Elmore and Sands 1949), optical 


resolution, target chamber vacuum and an approximate wavelength scale 
are indicated. 


ee ee 8 eee 
} Figures 1 and 3 are shown as plates. 
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Table 1. Emission lines from solid potassium 
(Wavelengths for M,,V are for the edge and extrapolated low energy limit) 


Peat Penieteee | Hales Inst. | Peak intensity 
Line a ye a e ae width | counts/sec/ma 
A) (ev) (ev) ane 
(ev) at 4 kv 
OTT9, 18-28, 
+0°5 +0-01 
M,,\ ae 0-025 500 
744 16-66, * 
+1 +0-02 
L,M, 47-63, 260-2, 1-96 
(Ll) 40-15 +058 40-1 eS Sw 
L,M 47-22, 262-5, 1-96 
( 1-4 1050 


Additional information: 


Ey, — Ey = 2°25 + 0-3 ev. 
Edge width for M,,V=0-19+0-01 ev. 


§ 4. DIscussION 
4.1. The M Valence Spectrum 


The M band emitted when conduction (4s) electrons make transitions 
into ionized M,,(3p) levels is shown in fig. 1 and has been observed 
previously by Kingston (1951) and Fisher et al. (1958). 

The b.c.c. alkali metal, potassium, has one valence electron per atom 
and is therefore expected to have a half-filled first Brillouin zone. In 
addition, since the first zone for an ideal b.c.c. lattice may contain 
1-48 electrons per atom before the Fermi surface departs from a sphere, 
and since it has been calculated that for the real metal potassium the 
relative effective mass of the conduction electrons is 0-94 (Kittel 1956), 
it is anticipated that the conduction electrons in potassium will behave 
very much like free electrons. Hence the density of states function N(£), 
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is expected to vary parabolically with # as measured from the bottom 
of the band, to show a sharp cut-off at the Fermi surface and to show 
no evidence of zone overlap in the occupied states. 

This behaviour is observed in the M emission curve of fig. 1 and in fig. 2 
where the measured ordinates on the emission curve, which are recorded . 
on a wavelength scale, are multiplied by the factor A° and replotted on an 
energy scale to give the f(H)N(£) curve (Tomboulian 1957). The broken 
curve in fig. 2 represents a parabolic (or H1?) extrapolation of the low 
energy portion of the curve in order to ‘remove’ the tail (Skinner 1940). 
The extrapolated bandwidth so obtained is 1-62, + 0-04 ev which is smaller 
than the value 1-9+0-2ev given by Kingston (1951). The extrapolation 


Fig. 2 
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Plot of f (#)N(£) vs. # deduced from the potassium M emission band. 


shown in fig. 4 of Kingston’s paper is hard to justify, and gives an excessive 
value for the bandwidth though the emission curve given is in close 
agreement with the present finding. Both these experimental values 
are less than the Sommerfeld ‘free electron’ width of 2-lev. 

The observed width of the emission edge is 0-19 + 0-Olev from 5% to 
95% of the maximum intensity with the instrumental width 0-024ev. 
It would appear, therefore, that about 0-17 ev of the edge width must be 
accounted for by temperature effects and by the width of the core level 
(M, or M;) in emission. 

There is no evidence of the edge doubling which is seen in the L spectra 
of the metals magnesium and aluminium, and which arises from the doublet 
nature of the 2p levels. The absence of edge doubling may result from 
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either the M,, M,; separation being small (less than a few hundredths of 
an electron volt) or the M, band intensity being very low. Since an M,, M, 
separation of at least some few tenths of an electron volt is expected by 
extrapolation from elements of higher atomic number, it must be concluded 
that the M, band intensity is very low (<1% of the M, band intensity). 
This conclusion is not entirely consistent with the results obtained by 
Skinner et al. (1954) for the M,, M, separations and relative intensities for 
the iron transition series, but a similar near absence of the L, spectrum 
was observed in sodium by Skinner (1940). 

The 0-17ev breadth of the emission edge must then be explained partly 
by the width of the M, level and partly by temperature effects on the 
Fermi edge whose width from 5% to 95% of the maximum is equal to 6k7’. 
Since 6k7' has values of 0-15 ev and 0-17 ev at temperatures of 20°c (cooling 
water temperature) and 56°c, it is seen that the maximum possible value 
of the M, width is 0-02ev, or alternatively that the maximum possible 
temperature of the emitting atoms was 56°c. Either conclusion is quite 
consistent with the experimental conditions. 


4.2. The Diagram Lines and the L Valence Spectrum 


The various lines, other than the M valence band spectrum, are illustrated 
in fig. 3. The identifications are indicated and were made using the energy 
level values of Hill e¢ al. (1952) and Fine and Hendee (1956) and the 
wavelengths of the Ll, Ly doublet given by Tyrén (1937). The peak 
wavelengths, widths, etc., are listed in table 1, as well as data for the two 
broad unidentified lines which for convenience are designated X and Y. 


4.2.1. The 474 doublet 


The doublet at 47 A emitted when electrons make transitions from the 
M, (3s) level to the two L,3 (3p) levels, was unresolved in first order but 
partially resolved in fourth order. Separation of the observed fourth-order 
spectrum into its two components was effected by assuming two nearly 
Gaussian profiles of equal widths with the intensity ratio of 2-2:1. 
Though this differs from the ideal ratio 2:1 it is consistent with the rapid 
fall of grating reflectivity in this region near cut-off. The absolute 
wavelength of either component is determined to 0-15 A from fourth-order 
spectra, while the separation is measurable to within 0-05 A from the same 
spectrum. 


4.2.2. The 2844 line, LL, 


Emission lines resulting from transitions within the L shell have been 
observed by Skinner (1940) and identified by Tomboulian and Cady (1941). 
The line identified as L,L, in the present case is very weak (fig. 3), but 
it appeared on all records and its energy may be computed fairly accurately. 
Little can be said about the width except that it is less than 0-75ev. 
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The line is identified as L,L, rather than L,L, or L,L3 since the L,L, 
intensity is expected to be at least twice as great as that of L,L,, and ae 
is no sign of another line 154 distant as would be anticipated from the 
proposed L,, L; separation of 2-25 ev. F 

The analogous line M,M,, expected at about 6204, was not observed. 


Fig. 4 
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Conventional energy level diagram for solid potassium constructed from 
observed emission wavelengths. 


4.2.3. The Ly valence spectrum 


A faint line at 414A is identified as the L,; valence spectrum (fig. 3). 
The intensity is very low on account of the proximity of the grating 
reflectivity cut-off but the line was seen on all records. The peak 
wavelength was computed by measuring the separation from the Ll, Ly 
doublet. This procedure gave a probable error of +0-38A instead of 
+074 if the peak wavelength had been obtained directly from the 
calibration table. Since the instrumental width in this region is 7:2ev 
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a value is not quoted for the band width. However, the ‘line’ will be 
several electron volts broad since the 1-62 ev spread of the valence electrons 
and the 2-25ev separation of the L, and L, levels both contribute to its 
width. 
4.3. Unidentified Lines 

There remain two lines, X and Y, which will not fit in with the proposed 
energy level scheme of fig. 4. In addition, their widths are excessive and 
not compatible with the approximate level widths of fig. 4. 


4.3.1. The line X at 248 eV 


It is not feasible to ascribe the line X to emission from contamination 
of the target since its intensity was high, its intensity relative to Ll did 
not change with time and there was no trace of carbon K radiation or 
visible deterioration of the target after bombardment for half an hour 
or more. The line has also been observed by Jacob (1956, private 
communication). The line is not a satellite with an initial state of the type 
KL,, KLs, etc., since it is still observed at the same intensity relative to Ll 
when the target voltage is reduced to 2-5kv. The ionization energy of 
the K level is about 3:6kev. Furthermore there is no obvious parent line 
for a satellite with an initial state of the type L,L,, L,L,, ete. 

If it is assumed that one of the states involved in the transition giving 
rise to the line X is L, or L;, then the other must be a band of excitation (?) 
states, several electron volts wide and lying some 10 ev below the M, level. 
Such excitation states have been observed by Parratt (1959) for the 
valence levels of manganese. 


4.3.2. The line Y at 120eV 


As with the line X, the line Y cannot be fitted to the proposed energy 
level scheme. However, its intensity is very low ( ~ 1/20 that of the line X) 
and it could conceivably have been emission from a contaminant though 
such a possibility seems remote. It is too broad to be the result of a 
transition between two sharp atomic levels though it may be an unresolved 
doublet. Once again it is hard to identify it as a satellite line for lack of 
a suitable parent potassium line., 


§ 5. Toe ENercy Leven Diacram 


The conventional energy level diagram for solid potassium, constructed 
entirely from the present data, is shown in fig. 4. The transitions observed 
are indicated, as well as the energies and widths computed from them for 
the valence, and the three L and three M levels. In constructing the 
diagram the peak of the L,;V spectrum was assumed to be the peak of the 
L,V component which is expected to be at least twice as intense as L,V. 


§ 6. L, Levets or Sopium, MAGNESIUM AND ALUMINIUM 


The transition L,L, in the metals sodium, magnesium and aluminium 
has been observed by Skinner (1940) and identified by Tomboulian and 
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Cady (1941). The value of the L, energy level relative to the highest 
occupied valence level may be calculated by summing the energy of the 
L,L, transition, and the transition energy corresponding to the L, emission 
edge. Tomboulian and Cady (1941) have made this calculation using 
Skinner’s (1940) data. 

The L,L, transition has been observed in the present case for these 
three metals. However, while the wavelengths for sodium (379-6 + 1-04) 
and magnesium (317-8 + 0-5 A) are in fair accord with previous results, that 
for aluminium (274-5 + 0-5 A) is not. A reason for the discrepancy is not 
forthcoming, though Skinner’s (1940) wavelengths are quoted only as 
rough values in a footnote. The relevant data, and the calculation for 
the L, level energy, using Skinner’s (1940) wavelengths for the L, emission 
edges, are given in table 2. 


Table 2. Evaluation of the L, energy level in Na, Mg and Al 


Na Mg Al 


AL, edge 405-45 + 0-05 250-45 + 0-05 170-56 + 0-05 
(Wavelengths from Skinner 1940) 


Ey, edge (ev) 30-577 + 0-004 49-500 + 0-010 


AL, L, 379-6+1-0 317-8 +0-5 274-54 0-5 


(Wavelengths from present investigation) 
(375) (317-4) (290) 


co 32-66 + 0-09 39-01 + 0-06 45:16 + 0-08 
(33-06) (39-06) (42-67) 


L, energy 63-24 + 0-09 88-51 + 0-07 117-85 +0-10 
(ev) (63-64) (88-56) (115-34) 


The values in parentheses, without errors, are from Skinner’s (1940) wavelengths 
for the peak of the L,L; lines. 
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The Influence of Solutes on Self-diffusion in Metalst 
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ABSTRACT 


A statistical theory of the effect of solutes on solvent self-diffusion in dilute: 
alloys is given. It is based on the vacancy model used previously by 
Hoffman, Turnbull and Hart and includes a qualitative consideration of 
correlation effects. Two limiting cases are considered: (i) in which the vacancy 
jump frequencies are not greatly changed by the solute atoms (‘ weak 
perturbation ’ limit), (ii) in which the vacancies attach themselves strongly 
to solute atoms and where the jump frequencies may be greatly different 
(‘ tight-binding’ limit). In both cases it is shown how experimental data. 
on solute enhanced diffusion of solvent atoms combined with knowledge of 
solute tracer diffusion rates may be used to deduce ‘ correlation factors ’ for 
solute diffusion. Analysis of existing data on dilute Ag alloys yields similar 
correlation factors whichever of the two limiting equations we use. The 
results are not consistent with the absence of isotope effects in Cd and In 
diffusion in Ag: this may indicate that the model is inadequate for the 
discussion of solute enhanced diffusion. Further experimental work, 
particularly on isotope effects, is called for. 


§ 1. InTRODUCTION 


A number of authors (Seitz 1950, Overhauser 1953, Hoffman ef al. 1955, 
Hart et al. 1955, Hoffman 1958, Lazarus 1955, Lomer 1958, Reiss 1959) have 
treated the problem of the influence of small concentrations of solutes on 
metallic self-diffusion. These discussions have been based on a vacancy 
model of diffusion and have assumed the frequency of solute-vacancy 
exchange to be different from the rate of solvent-vacancy exchange in the 
pure solvent. The frequencies of solvent-vacancy exchange in the 
immediate vicinity of a solute atom also are assumed to be different from 
their values in the pure solvent. Longer range effects have sometimes been 
included (Overhauser 1953). In the present note we shall concern our- 
selves with the model which treats the effects of the solute atom on vacancy 
jump frequency as localized. This is the model which most authors have 
used, but their results are not always consistent. In particular, various 
results have been presented for the concentrations of vacancies in.a dilute 
alloy. In the next section we present the statistical thermodynamic argu- 
ment by which the correct results for the local model may be derived and in 
the following sections we consider the application to self-diffusion in a very 
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dilute alloy. Then in §3 we show how an analysis of the effects of solutes 
on solvent self-diffusion may be used to make inferences about solute 
diffusion. In particular we are able to derive values for the correlation 
factor for solute diffusion (LeClaire and Lidiard 1956, Manning 1959). 
Knowledge of these factors is necessary to theories of the elementary jump 
processes. The inferred values for a number of solutes in Ag are given in 
the table. Those for Cd and In are not consistent with the absence of 
isotope effects in the diffusion of these elements (Schoen 1958 a, b). This 
may indicate the inadequacy of the present model for the description of the 
effect of solutes on solvent self-diffusion ; as the isotope effect experiments 
in principle yield more direct information about correlation factors their 
extension to other solutes would be valuable. 


§ 2. Vacancy CONCENTRATIONS IN A DILUTE ALLOY 


We assume the alloy to be composed of N solvent atoms and JN; solute 
atoms randomly distributed. It will be supposed that the energy of forma- 
tion of a vacancy at a site which is a nearest neighbour of a solute atom is 
different from the energy of formation at more distant sites, but we neglect 
any interaction between solute atoms and vacancies which are not nearest 
neighbours (this is in general agreement with the electronic theory of these 
interactions). We therefore specify the state of the system by giving the 
number of nearest-neighbour solute-vacancy pairs and the number, ny, of 
‘free’ vacancies. Let the former be JN,p, so that a fraction p of the im- 
purity atoms are associated with vacancies. 

If the Gibbs free energy of vacancy formation is gy for a free vacancy and 
(gv + Ag) for an associated vacancy, i.e. — Ag is the energy of solute-vacancy 
association, the free energy of the crystal containing ny free vacancies and 
Np associated vacancies includes a term 


nvgy + Nip (gv + Ag). te oh ete See 


We must also calculate the configurational entropy associated with the 
various arrangements of defects all of energy (2.1), To do this we note 
first that the total number of lattice sites is 


N,=N+N(1—p)+2N,p+ny=N+N,4 Nip t+ ny. . . (2.2) 
The number of ways of arranging the N,p solute-vacancy pairs among 
these J, sites is 
(arp P=) 


(N,p)! I (N+ N;+ Nip + ny— 2s). a er (acca) 


Here z is the number of nearest neighbours, i.e., the number of distinct 
orientations of a solute-vacancy pair. Expression (2.3) neglects small 
hindrances to the rotation of one pair caused by the proximity of others; 
these are negligible at the low concentration limit of interest here. After 
we have arranged the solute-vacancy pairs in one of the above (2.3) con- 
figurations there remain N+ N,—N,p+ny sites. The ny free vacancies 
N,(1 —p) free solute atoms and N solvent atoms are to be arranged on foe 
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sites subject to the condition that no vacancies and impurities be nearest 
neighbours of one another—since they would then constitute a solute- 
vacancy pair. First, the number of ways of putting the N,(1—p) free 
solute atoms onto the V + V,(1—p)+ny available sites is 


(N+N(1—p)+ny)! 

[V(1—p)]!(N+ny!)" 
There then remain NV + NV,(1—p) + ny less (2+ 1)N,(1—-p) sites on which to 
arrange the n, vacancies. This can be done in 


(N=2Ni(1—p)+ nv)! WA aed Pais 5) 

[N—zN,1—p)]!ny! 
ways. The product of the expressions (2.3), (2.4) and (2.5) gives the total 
number of configurations of energy (2.1). Actually in these configurations 
we have included some, e.g. a solute atom in the vicinity of a solute-vacancy 
pair, to which strictly we should assign an energy different from (2.1), but 
these configurations are so few that no significant error is introduced into 
the calculation. From expressions (2.1)-(2.5) we can immediately write 


down the free energy of the crystal: relative to a crystal without vacancies 
it is 


(2.4) 


Np) 


G=—kT In nf oe ae Tl (N +N, +N,p +ny— 2s) 


(WEN (1 —p)+ 1) 
[Vi(1—p)]!(NV+ny)! 
- (N —zN,(1—p)+nvy) ! 
[N —zN,(1—p)]|!nv! 
The average therodynamic values of ny and p are obtained in the usual 
way by setting 


} mg + Nips + Ag) . (2.6) 


0G 0G 

pete) eet 

on Op 
It is not difficult to show that these conditions lead to the following results, 
valid to first order in the concentration of solute atoms ; 


i eee any: te ay tere 

N.—Np—2zNiAl—p) exp ( je) a 
and 

= _ 9vt Ag 28 

N,p=Neexp( RT ). PS) Gee ete (2.8) 


Hence the average number of free vacancies ny is the same as in a pure 
solvent with only N,—N,p—zN,(1—p) sites. The terms subtracted from 
the actual total number of sites V's are the number of associated vacancy 
sites (N,p) and the zN,(1—p) precluded sites around the unassociated 
impurities. The degree of association p, by (2.8), will in general be small 
compared to unity and (2.7) may then be re-expressed as 


ay =(l—2c) exp(—gv/RT) . . . . «=~ (2.9) 
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where %y=ny/N, and c=N,/N, are molar fractions of free vacancies and 
solute atoms. Equation (2.9) is the expression used by Lomer (1958): 
it has also been obtained by LeClaire (unpublished) who calculated the 
entropy term by writing down the number of ways that the N,p pairs 
could be arranged on the N 2/2 bonds of the lattice. In the following 
sections we shall use (2.7) as well as (2.8). 


§ 3. APPLICATION TO SELF-DIFFUSION IN DILUTE ALLOYS 
oF FACE-CENTRED CUBIC STRUCTURE 


The elementary random walk theory of diffusion in an isotropic lattice 
gives the following expression for the diffusion coefficient (see e.g. 
Chandrasekhar 1943), 


DoD. ea ee eee 


in which I is the average number of jumps made by a diffusing atom in unit 
time and s is the magnitude of the jump distance. A more complete 
analysis of vacancy diffusion includes a consideration of the effects of 
correlation between successive vacancy jumpst (Bardeen and Herring 
1952). Inthe present connection we do not aim to include these rigorously ; 
we proceed first on the assumption that (3.1) is valid and indicate possible 
ways in which correlation effects may enter. We confine ourselves to a 
face-centred cubic lattice and begin by calculating [. 

In the present model of a dilute alloy there are several distinct vacancy 
jumps; these and their associated frequencies of occurance must be con- 
sidered separately to arrive at a correct expression for I for the solvent 
atoms. We use the following notation taken from earlier papers on solute 
diffusion: an associated vacancy exchanges places with its solute atom w, 
times/sec, it exchanges places with a given one of the four sites which are 
nearest neighbours to both the vacancy and the solute w, times/sec and it 
jumps to a given one of the remaining seven of its nearest neighbour sites 
(i.e.dissociates) k, times/sec. The corresponding association frequency from 
a next nearest neighbour site of the solute to a nearest neighbour site is ky. 
The frequency with which an unassociated vacancy jumps to another 
unassociated site is called wo. 

With this notation and the result (2.8) we can write down the number of 
solvent-vacancy exchanges/sec made by vacancies initially in positions 
nearest -neighbour to a solute atom ; 


kT 


We also want the number of association jumps which occur. Now, 
surrounding any solute atom there are 42 sites from which a nearest 
neighbour position can be reached in one jump; altogether there are 84 
a 


12 Nyexp(— 25572) (aun, + Th). a ichdcin Poin (BIE) 


} This is the original paper on the topic : i i 
ee Ape per on the topic ; for a brief review of subsequent 
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distinct jumps, each of frequency k,. We want the number of impurity - 
vacancy pairs at these separations. Now the probability that any lattice 
site is occupied by a vacancy is ny/N, if no solute is in the surrounding 
nearest neighbour shell. Also the probability that any site is occupied by a 
solute atom is N,(1—p)/N,, if no vacancy is a nearest neighbour of that 
site. Hence the total number of solute-vacancy pairs at the separation 
from which an association jump can be made is 


42 Ni(1—p) ny/N,. 
The total number of association jumps occurring per sec is likewise 
84k, Nj(1—p) ny/N,. ee Mrer ent eee GID) 
The remaining 420 jumps possible from these 42 sites around each solute 


atom are Ww») jumps and we shall assume that these occur with the same 
frequency as in the pure solvent. Their contribution is 


226 weN (pln NeeO Ce oe oe 2h S(94) 


So far we have accounted for 42N,(1—p)ny/N, of the total of ny free 
vacancies. The remainder, 


Ny(1—42N,(1—p)/N,), 
execute a total of 


L2wyty(1—42N{l—p)/N,) . ©. . « . (3.5) 
jumps per sec. We may note that inclusion of terms (3.4) and (3.5) is only 
possible for solute concentrations less than about 2%; otherwise over- 
lapping of the spheres of influence of different solute atoms becomes 
important. The total number of solvent atom jumps per sec is given by 
the sum of terms (3.2)—(3.5) with ny given by (2.7), and I’, the average 
solvent atom jump frequency, is obtained by dividing this sum by the 
number of solvent atoms present. To the first order in c the result is 


T= 12w, exp ( —guikD)| 1—18c+cexp (— agin) | (3.6) 
0 
in which we have used the result that k,/k, =exp (— Ag/RT). 

When now we wish to consider the diffusion coefficient of the solvent 
atoms it is necessary to include correlation effects. The simplest situation 
to handle is that where the solute atoms may be regarded as constituting a 
‘weak perturbation’; then none of the jump frequencies will be changed 
radically and correlation effects will be closely the same as in the pure 
solvent. Under these conditions we can go immediately from (3.6) to the 
equation 


4 14k 
Ds(c) = D35(0) [oe 18¢+cexp (— Ag/kT) Gree ‘) | (3.7) 
0 
in which D,(c) is the diffusion coefficient of the solvent in an alloy containing 
a concentration c of solute. 
Equation (3.7) may be a reasonable starting point. We have at present 
little knowledge of the relative magnitudes of w,, We, ky and wo, but a 
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‘weak perturbation’, assumption seems a better approximation than a 
model of tightly bound solute-vacancy pairs (§ 6), so far as one is able to 
judge directly from the relative diffusion rates of solutes and of the solvent. 
Other possibilities will be noted later. 
Equation (3.7) may be written 
4(w, + 7k,/2) Ag 

D,(c)=D,(0) (1+bc), b= —18+ Serger aie (3.8) 

Now, this same model of the alloy yields the expressions 


_ 2s®w, (w, + Tk, /2) ( _ gvt =!) 3.9 

is (w,+w,+ Tk,/2) we kT aera kh, 
for the solute (‘impurity ’) diffusion coefficient and 
qv 

D,(0)=25%(9/11) exp (— £) : Rae, et): 


the 9/11 being an approximation to the correlation factor, fy (consistent 
with (3.9)) for solvent self-diffusion in the pure metal. The correlation 
factor for solute diffusion into pure solvent is the factory 
f,= (wy + Thy]2)|(wy + wy t+ Thy/2) . . . . (3.11) 
n (3.9) (LeClaire and Lidiard 1956); this is a quantity of interest to the 
electronic theory of solute diffusion since only if it is near to unity is the 
solute diffusion determined by w, alone, as often assumed. It is therefore 
interesting to note that eqns. (3.8), (3.9) and (3.10) allow us to write 
4f D,; 
=1— : |, 41 Rey Se Se a 
fm Get (Baa) soe 
We may note that this is consistent with the description of a weak perturba- 
tion, for as D,/D,(0)>1 and 6+0, then f;>f,=9/11, as it must do. 
The expression for 6 corresponding to (3.12) is 
36 D, 
wee: (575) By Ahem Caen, rds 
A qualitative consideration of correlation effects in other than the ‘ weak 
perturbation’ limit indicates that (3.12) changes mainly in that the 18 is 


reduced. For a number of solutes b is several times 18 and the results are 
not greatly affected. 


b6=—18+ 


§ 4. ExpurmMenTAL RESULTS FOR Ag ALLOYS 


Experimental results for a number of Ag alloys are available. D5(c) is 
often well represented by an exponential dependence exp (bc) rather than 
by the linear form (3.7); however, we are only interested in the limit of 
c<2% and it is reasonable to assume that 6 in the exponential form has 
the same significance as here. Results for f, calculated from (3.12) for 
1000°K are given in the table, column W.P. The use of Manning’s 


EEE naseeeere ee ee 
} A more accurate expression for f; is available in th i 
RnR RRA i e special case when 
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expression in place of (3.11) and the consistent assumption k, = Wy decreases, 
the values of f; in the table only very slightly. 


Values of 6 and D; D(0) for various solutes in Ag at 1000°K (as interpolated 
from experimental results) and the values of the solute atom correlation 
factor f; as deduced in weak perturbation and tight-binding 


approximations 
Solute b D jj D,(0) gale Bae tip AAs 
Pd — 8-2? 0-123» 0-96 eS 
Cu 10-2° 1-254 0-86 Fa 
Cd 6-9° 367° 0-52 0-46 
In AR hs 530° 0-43 0-48 
fi! 36-46 8-39f 0-50 0-72 
Ge 32-7° 9-51! 0-39 0-61 
Pb 86-0° EAS 0-65 0-89 
Sb 51-78 7-658 0-64 0-86 


Sources of experimental data in this table : 
(a) N. H. Nachtrieb, quoted by H. Reiss in Phys. Rev., 118, 1445 (1959). 
(b) N. H. Nachtrieb, J. Petit and J. Wehrenberg, J. chem. Phys., 26, 106 

(1957). 
(c) R. E. Hoffman, D. Turnbull and E. W. Hart, Acta Met., 3, 417 (1958). 
(d) A. Sawatsky and F. E. Jaumot ; quoted by R. E. Hoffman, Acta Met., 
6, 95 (1958). 

e) A. H. Schoen, Ph.D. Thesis, University of Hlinois (1958). 

f) R. E. Hoffman, Acta Met., 6, 95 (1958). 

g) E. Sonder, Phys. Rev., 100, 1662 (1955). 

h) E. Sonder, L. M. Slifkin and C. T. Tomizuka, Phys. Rev., 98, 970 (1954). 


For the calculation of D;/D,(0) we have used 
44-27 k cals 
D,(0) =0-442 exp (- ce cm?/sec 
as given by Schoen’s analysis of the data of Slifkin, Sonder and Tomizuka. 


§ 5. Discussion oF RESULTS IN § 4 


The correlation factor for solute diffusion in pure solvent, f,, must by its 
definition lie between 0 and 1. It is encouraging that the values in the 
table all do lie within this range. It is encouraging also that the value for 
Cu in Ag, which we might expect to be a closely ideal system, is close to the 
value for Ag self-diffusion (i.e. 0-82 in the present approximation). There 
are indications that this is true also for Au in Ag (Slifkin, unpublished data). 
It is also worth drawing attention to the fact that none of the /; values is 
small, and, although this is not a sensitive test of our ‘weak perturbation ’ 
assumption, it is nevertheless consistent with it. 

On the other hand, it must be pointed out that there exists the quite 
independent, and apparently direct, evidence of the isotope effect experi- 
ments for the diffusion of Cd and In in Ag which shows that the correlation 
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factors for these solutes are small, i.e. that their diffusion is controlled by 
(w,+7k,/2) and not by w, (Schoen 1958a, Tharmalingam and Lidiard 
1959). If these results are confirmed and extended to other elements we 
must infer from the table that a considerable part of b, i.e. of the enhance- 
ment of self-diffusion, is due to longer range effects than we have considered. 
Thus, while it appears from Schoen’s work that the Overhauser type of 
theory cannot account for all the enhanced Ag diffusion caused by Cd 
and In solutes (Schoen 1958b) it may account for part of it and so give 
lower f; values than we have obtained here. 


§ 6. Comparison with Work oF HorrMaAN, TURNBULL AND Hart 
AND WITH ‘ TIGHT-BINDING ’ MODEL 


These authors (HTH) made the first quantitative application of the 
model described in this note to the problem of the enhancement of metallic 
self-diffusion by solutes. Not surprisingly, their result is similar in general 
form (see Hart et al. 1957). One obvious difference, though, is the occur- 
rence of the number of nearest neighbours, 12, in place of the 18 in our 
eqn. (3.13). Thus their equation in our notation and for the f.c.c. lattice 
reads 


S79 or 
b= 12+8, (575), 
in which a, is their quantity, “the number of effective solvent diffusion 
jumps per solvent atom while in the first co-ordination shell of the solute ”’. 
The 12 in eqn. (6.1) arises in the HTH treatment through the use of eqn. (2.9) 
and the assumption that this gives the fraction of time that a given solvent 
atom diffuses with diffusion coefficient D,(0). The two treatments come 
closest to formal similarity if we take k,=wp, i.e. the assumption of very 
short range of interaction which is sometimes made. Terms (3.3) and (3.4) 
may then be taken together and the constant term in 6 then becomes — 11 
instead of — 18 as in (3.8) and (3.13). The remaining difference of 1 is due 
to an apparent neglect by HTH of the fact, included in (3.6), that the molar 
fraction of solvent is 1—c and not 1. Actually such a comparison is rather 
forced since the HTH treatment, although it can be said to include the 
association jumps (if k,=wp»), does not explicitly include the dissociation 
jumps (‘,). In fact if k,=wy, then k, exp (— Ag/RT) also equals wy and in 
place of (3.8) it would be consistent to write 


4w A 
i Sih Der, = g 
or ae exp ( a) Sete uA An er OS) 


The real point of alternatives such as these must however lie in the corre- 
lation effects. The correlations between subsequent jumps of a solvent 
atom and a prior jump are unlikely to be very different from those in the 
pure lattice if the prior jump is either a k, or a k, jump, particularly if 
k,=Wwot. On the other hand, if the prior jump is a w,-jump then the 
a 


(6.1) 


} Dr. J. R. Manning in a private communicati is j . 
bearers Pp on reports that this is confirmed 
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correlation effects may be different from those in the pure lattice especially 
when w,<w,. If we introduce a different correlation factor for the 
w,-jumps we must replace (6.2) by 


b= 44 Troxp(— of), (hs 20G) ite x inc) tes a 053) 
Wolo 

in which f, is the correlation factor for the pure solvent metal (=0-78 
accurately); the second correlation factor f,= f,(w./w,) is for an initial 
jump w, and has been calculated by Compaan and Haven (1956) and, more 
approximately, by Reiss (1959). Equation (6.3) and the corresponding 


expression 
4f iD: 
b= —-—4 z : Aue ; 
| + =a (Davy): pd 
with 
f= 1|/(W1 +2), . . . . c . (6.5) 


may be regarded as reasonable approximations to the ‘tight-binding’ 
limit, i.e., the opposite limit to (3.8) and (3.13). Owing to present limited 
knowledge of the function f,(w,/w,) solutions of (6.4) and (6.5) must be 
obtained by graphical interpolation between the values given by Compaan 
and Havan (1956). This yields the f, values given in the table under the 
coulmn headed T.B. (the second decimal figure is uncertain). For Cd 
and In the figures are little altered from those obtained in the weak per- 
turbation limit; for the elements Tl-Sb they are increased by about 0-2. 
These changes come primarily from the incorporation into the equations of 
J,; the values of which are several times less than the 9/11 assumed in the 
weak perturbation treatment. This means that in order to give the same 
observed effects it is necessary to have larger ratios of w/w, than are 
obtained if we use eqn. (3.8) (cf. the product w,f, which occurs in (6.3)). 
This in turn leads to larger values for /;. 

In a recent private communication to the author Dr. J. R. Manning 
reports a more extended discussion of the model with k,=w,. He does not 
neglect k, in comparison with w, and w, as we have done. In addition 
he uses a more accurate expression for f, (Manning 1959) and develops 
expressions for f, as function of w,/k, and w,/k,. Experimental data on 6} 
and D,/D (0) is then used to infer w,/k, and w,/k, and hencef;. The inferred 
values of w,/k, and w,/k, are large compared to unity and the f; values agree 
well with those in the table, column T.B. These results may mean that 
more reliance should be placed on the T.B. values. However the inference 
of separate values for w,/k, and w,/k, is only possible on account of the 
assumption that k,=w,; the sensitivity of the analysis to this assumption 
is not known. 


§ 7. CONCLUSION 


The table shows that the inferred correlation factors are not highly 
sensitive to the alternative mathematical assumptions which we have made 
about relative jump frequencies entering into this problem. In particular 
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the values for Cd and In seem to be firm. These values are however in 
conflict with the much lower values demanded by Schoen’s isotope effect 
experiments. The inference may be that the present model is too simplified 
to deal with the problem of solute enhanced diffusion and that effects of 
longer range must be included. However, we emphasise that the analysis 
presented here is perfectly consistent considered on its own, and that any 
such inference rests entirely on the isotope effect results for Cd and In in 
Ag. Jt would clearly be of interest to have these measurements extended 
to other solutes. 
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ABSTRACT 


The mathematical methods used in a calculation of nuclear interactions 
by the optical model are described, and the most efficient way of programming 
the calculation for an electronic computer investigated. 


§ 1. GENERAL FORMULATION 


WHEN a nucleon or group of nucleons interacts with a nucleus, the 
behaviour of the system is determined by all the possible interactions 
between the nucleons present. This is very difficult to handle mathe- 
matically, so the optical model replaces all these separate interactions by 
a two-body interaction between the incident particle and the ee nucleus 
of the form (Bjorklund and Fernbach 1958) 

V(r) = Ve(r) + Uf(r) +tWog(r) + (Us +iWs)h(r) lio eel) 
where yr is the separation of the two particles, and U, W, Us and Ws 
are the real and imaginary parts of the central and spin-dependent 
potentials respectively. V.¢ is the Coulomb potential, f(r), g(r) and h(r) 
are functions giving the radial variations of the potentials (cf. §4). This 
is equivalent to replacing the target nucleus by a potential well having 
both refracting and absorbing parts. The refracting potential U elastically 
scatters the incident particles, and the absorbing potential W takes 
account of all the inelastic processes, including inelastic scattering and 
absorption. The spin-dependent term given above is appropriate for 
the scattering of spin 4 particles by spin 0 nuclei; for the scattering of 
particles of zero spin there is of course no spin-dependent term, while for 
the scattering of particles of higher spins more complex forms are required. 
In all optical model calculations it is usual to take the spin of the target 
nucleus to be zero. 

To calculate the observable cross sections for two particles interacting 
through the above potential, Schrodinger’s equation is applied to the 
whole system 

2 
5, +H VyF=0 nage cc) 
where £ is the centre-of-mass energy. 
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If the total wave function ¥’ is expanded in spherical harmonies 
4 => THO) S CHT Sacer an 
jlm le 
we obtain a series of equations for each radial wave function. The spin 
of the incident particle can couple in two ways to the orbital angular 
momentum 1 to give the total angular momentum j=/+ 3, and the 
eigenvalues of 1.6 corresponding to these two spin orientations are l 
and —(l+1). The corresponding radial wave functions U;*(r) and 
U,-(r) satisfy the equations 
d?Ur*(p) [ Ve _ Uf(p)+iWg(p) 
tbo SS 


dp? E 
stiWs)h i+1 
2 WU, — )h(p) im” ( = | U,+(p) =0, ot a 
@UH(p) , fa Ve _ Uflp) +iW9(e) 
dp? E E 
gee teas SI ast Sa) oe = A Ur e)=0 (4b) 
P 


where p= kr, k being the wave number k= 1/(2u)/h where pw is the reduced 
mass of the incident particle, ~=A,/(A;+A,), A; and A, being the 
masses of the incident and target particles in a.m.u. These radial wave 
functions are zero at the origin and beyond the nuclear field tend 
asymptotically to the forms 


U,*(p) =F (p) + 1Gy(p) + S* LFi(p) -—1Gi(p)] - . - (5) 
where the S*+=exp(276,) are the (complex) elements of the scattering 
matrix, and F,(p) and G(p) are regular and irregular solutions of the eqn.(4) 


without the nuclear terms. These functions are taken to have the 
asymptotic forms 


F(p)——sin [p — y In 2p — la +04], Sew? paso Laas 
p> co 

G(p)——>-cos [p — y In 2p — 41a + ]. t+ EE 
p> oo 


The non spin-flip and spin-flip scattering amplitudes are then 


A(6) =fe(8) + a 2 {(1+ 1)S;+ +18,-— (21+ 1)} exp (2io,)P,(cos 6) (6a) 


Ne Be ‘ 
_ BO) = — > {S8;+—S}exp (2io,)Pp(cos#)  . . . . . . (6b) 
2k imo 


_ where fc(9) is the Coulomb scattering amplitude 


6 ; : ; 
fe(0)= — 57 Coser 5 OXp [2tog—2tyInsin#/2]  . . . (7) 
where ot ply Zine" 
ee” 


Z, and Z, being the charges of the incident and target nuclei, and o; is 
the Coulomb phase shift for the th partial wave. 
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The observable quantities may be expressed in terms of the scattering 
amplitudes. The differential cross section for elastic scattering is given by 


o(8) =|Al?+|Bl? ier g tet es aa, 

and the polarization by 
(6) = 2 Im AB* C 
aP+iBe" cfs Ve eon es!) 


where n=k, xk /|k, x k,l. 
The total absorption cross section is given by 


o4= Be > (E+ DU- |S") +1 |S-P)} ere (10) 


If the incident particles are uncharged, V¢, f¢(#) and a; are all zero, and 
the total elastic cross section is finite 


ee= Fe S48 HIS. (Y 
The total cross section is 
» co 

on= FD {(l+ (1+ Re S:*) +1(1— Re S,)}. ec. Oey 


In practice, partial waves with / greater than some value l,,,, do not 


contribute significantly to these summations and so they may be neglected. 


§ 2. MetHop OF SOLUTION 


The essential part of the calculation is the numerical integration of the 
radial wave eqns. (4) from the origin out to a distance p,, where the nuclear 
field is negligible. The wave function obtained by this integration of 
the internal equation must join smoothly to the external solution (5). 
Thus, if f,is the logarithmic derivative of the internal wave function at p,,, 


po OHA oles @—Oiel) an 
t  \ File) +1Gi(p) + SAL i(p) —2Gi(p)] J o= oy, 

Since f= is determined from the internal integration, and the F’s 
and G’s may be calculated as described below, this determines the S,*, 
and hence the cbservable quantities by the relations (6)—(12). 

The internal integration was initially coded for the Runge-Kutta-Gill 
method, since this is available as a subroutine on the computer used. 
Advantage can, however, be taken of the absence of first derivatives in 
(4) to use the Fox-Goodwin method 7, and this is up to ten times faster. 
In the Runge-Kutta method the value of the function and its derivative 
is required at one point beyond the origin and in the Fox—Goodwin 
method the values of the function at two points beyond the origin are 
required before the numerical integration can begin. If these points are 
taken close to the origin a series solution of (4) may be obtained for small p, 
and very few terms suffice to give accurate starting values of the wave 
function and its derivative. 

In the external region, the wave eqn. (4) reduces to 


oa +[ 12h SE | u=0. Repeal y hoo (TA4 
dp* pater <P 
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For y 4 0 its solutions are the regular and irregular Coulomb wave functions 
F(p) and G(p). The F, may be readily calculated by several methods 
among which the asymptotic series expansion, the Shepansky—Butler 
series and Miller’s method (Stegun and Abramowitz 1955) were investigated. 
Miller’s method was chosen for its simplicity and speed. This method 
sets F,, (p)=0, aF(p)=«, where / is an integer >1,,,, (-=60 was used 
here), « is a constant to be determined, and « a very small number. Using 
the recurrence relation 


[y2+ (t+ 1)? 7? ae Sey y 1 
Way eee Fyys(p) + ‘Ah F_y(p) = (20+ 1) | aon sm ~ |File) 
(15) 


all the «F;(p)’s down to af,(p) may be calculated. The constant « was 
found from the Wronskian 


Fq'(p)@o(p) — Fo(p)G@o (p) = 1 + Ny Seay 
where F’,'(p) was found from the recurrence relation 
y_ ttl ver Lay Ee 
En Shey Nps ey hy pry pL 
[St | Ale)—Fi(e)= 
and Go(p) and G,'(p) were calculated as described below. 


The irregular functions were evaluated using the expression (Froeberg 
1955) 


Fia(p) (17) 


Go(p) =s cos (— yn 2p +p +o 9)—tsin(—yIn2p+p+o,), 7 
G(p)=Se08(—yln2p+p+09)—Tein(—ylndptpto),f > 8 


where 
25 


25 25 25 
$m) 8h, whe tgs SS > ued ee 
n=0 n=0 n=0 n=0 
The number of terms in these sums was found by experience to be 
adequate, and the individual terms were calculated from the recurrence 
relations 


(nly 3c mee) 


Cee pM The yy 
me Qp(nt+1)™  ~2Qp(ne1) 
poe oy y?—n(n+ 1) 
ati oat So 
2p(n + 1) 2p(n +1) 
2n+1 ae (19 
gi fey eee Ost hye ) 
2p(n+ 1) 2p(n +1) p 
€ ' jt 
Eo (20+ ly ya n(n +1) ge tnt 


< ‘A: 
2p(n+1) ™  2p(n+1) ~" pp 
and the initial values 


B22 2 2 
$= 1, 4=+; S)=0, S\= B vs —2; 
2p 2p 2p 
y : 
t,=0, Pi ett LE ga Te ne ee 
2p p p» 2p 
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The summations satisfy the relation 
sT = St=1. Fete TRO BI ets LO) 

The convergence of these series improves as p increases, and to obtain 
accurate results a value of p substantially greater than p,, must be used. 
To ensure accuracy, p was first set to a suitable starting value, and the 
series summed. If the relation (20) was satisfied to one part in 104, the 
values were accepted; if not, p was increased until it was. 

The values of G(p;,) and G)'(p,;,) were obtained from these G,(p), 
Go (ep), where p>p,,, by integrating the eqn. (14) inwards numerically. 
The G(p;,) and G’(p,;,) were then found from (17) and (15). 

The relations (15) and (17) are satisfied by both the F’s and the G’s. 
They may only be used downwards in / for the F’s and upwards in / for 
the G's, when />p; otherwise instabilities develop and all accuracy is 
lost (cf. Hartree 1952). Likewise, eqn. (14) may only be integrated 
outwards for the F’s and inwards for the G’s in the unstable regions. 

The Coulomb phase shifts o,(y) were calculated from the asymptotic 
series expansion for large / (/=50 was used) 


of(y)=a(l+3)+ylnp—1] 
1 [ ‘ sinc sin3a sin 5a sin 7a sin 9a | (21) 


12 * 36082 126084 ~ 16808 11888" 


where 
a=tan[y/(l+1)], B=ly?+(l+1)}" 
and the recurrence relation 


Oy) = a(y)—tan™ [y/+1)]. - - - + (22) 
If y=0, the solutions of (14) are spherical Bessel and Neumann functions 
Fi(p)=pjilp), Gilp)=—pm(p). - - ~~ ~ (23) 


The F’s were calculated by Miller’s method and the G’s, using (19), with 
the initial values 
sin p 


UE ape , jrilp) = pe ae . > 
COS p cosp cosp 

1(p)= ——, LU a ager horace ot 

o(P) D 1\P 0 p 


The F’s and G’s and their derivatives were required to satisfy the 
Wronskian 
FY (p)Gi(p) — Fi(p)Gr (p) = 1 she etc sen (ae) 


for all / to two parts in 10°. 


§ 3. CHOICE OF PARAMETERS 


Four parameters determine the numerical part of the calculation: (a) 
the matching radius p,,, (b) the number of steps in the Coulomb integration 
Ne, (c) the number of steps in the nuclear integration N,, and (d) the 
number of partial waves (J,,,,+1). Within certain limits, the larger 


m 
these numbers are the more accurate the final result, but the longer the 


P.M. 4M 
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time of computation. It is therefore necessary to choose them to minimize 
the time consistent with the desired accuracy of the final results. 

Since the Coulomb integration is done only once, Ne can be chosen 
generously without substantially affecting the total time. This leaves 
py, Nn and l,,,, to be studied in detail. In general, such parameters 
can either be chosen once and for all, or calculated for each run by some 
empirical formula, or fixed automatically to satisfy some pre-determined 
accuracy criterion. The latter procedure is inherently the most satis- 
factory, but is only convenient for /,,,,. in the present work. The nuclear 
integrations were performed for increasing /-values, until 1—|ReS|| 
for the last partial wave became less than some small number e. Such 
an automatic procedure would be possible to fix p,, and Ny also, but it 
is not efficient as the nuclear integrations would have to be done again 
and again for increasing p,, and NV, until some pre-set condition was 
satisfied. It is easier to use an empirical formula for them. 

The best choice of p,, and NV, was investigated by calculating the 
scattering of 14Mev neutrons by aluminium for a range of values, and 
recording the accuracy with which sample cross-sections and polarizations 
were determined. For this investigation, /,,,. was chosen considerably 
larger than usual. Some typical results are shown in tables 1-3. Since 
most of the computation time is taken up by the numerical integration, 
and this is proportional to N,J/,,,,, it is apparent from the table that, 
adequate accuracy for most purposes can be obtained by choosing C=7 
where 


Pu =PRtCa . = . s * . = (25) 


where pp and « are the nuclear radius and diffuseness in p units, pp, =/R, 
a=ka. 

The best value of p,, chosen in this way is very insensitive to Z and A, 
so the above relation was written into the programme. The best values 
of N, depend on # and A and it was found adequate to have ten steps 
per unit of p. 

It would be possible to use varying step lengths for the numerical 
integration, the shorter steps where the wave function changes rapidly 
and the longer steps elsewhere. This would reduce the calculation time, 
but as the improvement is marginal it has not been done. A value of 
«=0-0001 was found to give accuracies at least as good as the best in 
tables 1-3. 

Combinations of the physical constants occur in three places in the 
calculation : 
the Coulomb constant 


s pL Le” 
kh? 


where £,,,, is in Mev. 
the wave number 


= 0-157454Z,Z,(A,/E,,)"2  . . . (26) 


k= 0-218728 (uiyq,)¥2 Scan eo Tr} 
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and the coefficient of the spin-orbit term (see § 4). 


TN ae 
(sa) = 0:011057 x 10-26 em? Be ig (28) 


Table 1. Percentage error in calculation of total cross section for 14 Mev 
neutrons on aluminium for various C and Ny, 


Number of integration steps Vy 
C 
15 20 25 30 
6 0-05 0-03 0-03 0-03 
% 0-06 0-02 0-02 0-02 
8 0-07 0-02 0-011 0-008 
9 0-11 0-02 0-008 0-003 


Table 2. As table 1 for differential cross section at 0° 


Number of integration steps Vy 


The programme was arranged to read the parameters of each calculation 
from a separate data tape. The parameters were A,, A;, Z>, Z,, the masses 
(in a.m.u.) and charges of the target and incident particles respectively, 
the energy £ of the incident particle in the laboratory system in MeV, the 
central and spin-orbit potentials (in Mev) and the parameters of their form 
factors (in 10-13 cm), the Coulomb radius r,, and the angular range required 
in the output. A second data tape provides the experimental data with 
‘which the calculation is to be compared. 


P.M. a 
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The output consists of all the input numbers, the range of the nuclear 
integrations, k, y, the S;*, the o,! and o4 (also of, og, of for Z,=0) 
do/dQ and P for the angular range specified, and A, and A> for each 
point of the experimental distribution, and their totals. 

The time in seconds for a calculation on the Ferranti Mercury computer 
is approximately 


T 20-5 H,,,, R? + 20 et ae ee 
where Ris the nuclear radius in fermis. This is halved if there is no 


spin-orbit potential. 


§ 4. RapiaL VARIATION OF THE COULOMB AND NUCLEAR POTENTIALS 


The Coulomb potential was taken to be that due to a charged sphere 
of radius R,=r,Al® 


Viz a. (3-7?/R.2) for r<R, 
~ ee for tS he a LT 
r 


It would be more accurate to use the potential due to a charged sphere 
with a diffuseness corresponding to that of the nuclear potential. This 
would substantially increase the computation time, without significantly 
affecting the results. The interaction is insensitive to the radial form of 
the potential and the form (30) is adequate for all practical purposes. 
This is shown by the comparison in tables 4 and 6 between two calculations 
with identical parameters except for r.. 

The nuclear potentials can be represented by a number of form factors 
all having a rapid fall-off beyond the nuclear surface. Those used in 
the present programme were the Saxon—Woods form 


(ie ge a ee 


1+exp (=) 
a 


where k=r)A"® is the nuclear radius and a the diffuseness of the nuclear 
surface, and the Gaussian form 


str) =exp| -(—*)"|. Se es 


The form factor for spin-orbit potential, by analogy with the Thomas 
term in atoms, was taken to be 


ho \* 1 G(r) 
bir) —'(—— | 6 = ge Ge 
”) (sa) ror i 3 

The factor (%/2Mc)*, the square of half the Compton wavelength of the 
proton, is introduced simply to keep A(r) dimensionless. Some authors 
use (//m,c)?, where m, is the pion mass. The corresponding potentials. 
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must then be divided by (2M/m,)?=180-59. The negative sign is 
introduced to keep h(r) positive. 

The programme was arranged to use the Saxon—Woods form for U, 
either the Saxon—Woods or the Gaussian for W, and the form given by (33) 
for the spin-orbit potentials. 

The choice of a definite p,, for matching the internal and the external 
wave function is equivalent to cutting off the potential at that radius. 
The problem is then solved exactly for the cut-off potential. Since little 
is known about the nuclear field at such distances from the nucleus, nothing 
is lost by this procedure. 


Table 4. Effect of altering Coulomb field. Calculation with same 
parameters as table 6 but with r,=1-4¥F (F=10-" cm) 


§ 5. PRoGRAMME CHECKS 


With a calculation of this complexity it is essential to check that no 
analytical or programming errors have been made. In the present case 
there are no simple limiting cases that could be used, and the whole 
calculation would take a prohibitively long time to check on a desk 
computer. The best check is to compare the results of some trial calcula- 
tions with that of similar runs with exactly the same initial parameters 
using another programme written independently for another computer. 
Through the courtesy of Professor M. A. Melkanoff such runs have been 
made, and the results are given in tables 5 and 6. The agreement is 
generally satisfactory, and the remaining small discrepancies may be 
accounted for by slightly different choices of the numerical parameters in 
the programmes, and to the accumulation of small rounding errors. The 
agreement is sufficiently good for comparison with experimental data. 


§ 6. COMPARISON WITH EXPERIMENTAL Data 


It is frequently necessary to compare the calculated angular distributions 
and polarizations with experimental data, and it is convenient to incorporate 
facilities for doing this into the programme. The experimental data forms 


4N2 
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part of the input, and when the rest of the calculation is complete, the 
values of 


ye Gee aoa apa5(209=P2)" ae 
ae do(8) on r= 2( 5Px(8) 
(“aa 


where the suffix H denotes experimental data, are calculated for each 
experimental point and included in the output. 


Table 5. Comparison of present programme with that of Melkanoff for 
17 Mev protons on gold, with U = —48 Mev, W = — 8Mev, rp= 13 F, 
re=1-3¥F, a=0-5F, A,;=1, A,=197. Both form factors are 
Saxon—Woods, and there is no spin-orbit potential. All cross 
sections are in millibarns per steradian, and are written in floating 
point form 


da/dQ da/dQ 
(Melkanoff) (present work) 


3-202, 4 3-208, 4 
1549, 3 1-549, 3 
2-951, 2 2-951, 2 
6-728, 1 7-624, 1 
1-460, 1 1-458, 1 
7-825, 0 7-807, 1 


984-7 mB 985-2 mB 


Table 6. Comparison of present programme with that of Melkanoff 
for 17Mev, protons on copper, with U=—46mev, W = —9mev, 
Us=—1444mev, Ws=361 Mev, 7r5=1-33F, re=1-33F, a=0-48F, 
A,=1, A,=64. (Im all these check calculations Melkanoff’s 
values of the constants in eqns. (26) and (27) were used. These 
are 0-158051 and 0-219438 respectively) 


| 9 do/dQ do/dQ P P 

| (Melkanoff) | (present work) | (Melkanoff) | (present work) 
20 2-807, 3 2-809, 3 —0-0590 — 0-0578 
40 1-819: 2 1-822) 2 0-0425 0-0426 
60 2-140, 1 2-136, 1 0-:2614 0-:2623 
80 1-563, 1 1-567, 1 —0:1213 —0-1168 
120 3°396, 0 3°389, 0 — 0:7876 —0-7811 
160 2-310, 0 2310.50 


0-6727 0-6646 


— 


oR 1054-7 mB 1055-0 mB 
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Experimental differential cross sections are often subject to a larger 
absolute than relative error. If the incident particle is charged, and the 
cross-sections are available to small angles, it is then desirable to normalize 
the experimental cross-sections to the theoretical in this region. A 
facility to do this for a specified number of experimental points was 
incorporated in the programme. 
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Giant Spin Density Fluctuations 


By J. C. PHituiest 
Royal Society Mond Laboratory, Cambridge 


[Received August 17, 1960] 


RECENTLY Overhauser (1960) has presented a one-dimensional model 
of interacting fermions which has the remarkable property that it 
exhibits in the Hartree-Fock approximation spin density fluctuations 
in the presence of an arbitrarily weak repulsive interaction. He has 
also argued that the qualitative features of his results can be obtained 
in two and three dimensions, and has thereby called into serious question 
the modern electron theory of metals. 

Overhauser’s spin-fluctuation state is based on energy gaps generated 
by exchange. The situation in one dimension is especially favourable 
for two reasons: 

1. Very little additional kinetic energy is required to generate the 
spin fluctuations. In the weak coupling limit (Overhauser’s no) the 
additional longitudinal kinetic energy vanishes exponentially. 

2. The exchange interaction is coherent in one dimension, since only 
one propagation vector q is used. In general the additional exchange 
energy AJ (we use Overhauser’s notation) contains a factor 6,,’. If the 
electrons are divided into m equivalent groups with different q’s, AJ is 
multiplied by 1/m. This has the effect of replacing n everywhere by mn. 

To see how these factors nullify Overhauser’s results in three dimen- 
sions, we consider the most favourable case, the cubic fluctuations 
sketched in Overhauser’s fig. 2. In contrast to the one-dimensional 
case the perturbed states here necessarily have an additional ‘transverse’ 
kinetic energy AU arising from the non-spherical occupied region. 
Again choosing the most favourable case, n=1, we find 


W ~2NE,(1—coth 2mn), 
~ —4NH,,exp(—12)=—0-024NH, . .. . (1) 
AU ~0-024NE,, Pie erage Lhe) 


so that in this case the non-magnetic, ferromagnetic and spin-fluctuation 
states all have nearly the sameenergy. Ifn <1, however, the ferromagnetic 
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state will be lowest (since its exchange energy is not reduced by the 
factor m-!) while if n>1 the non-magnetic state will be lowest. 


In normal metals 
n~4H, /(v)~10 as Seca 


with (v)=Ny/L~2ev, EH, ~5ev. Thus in normal metals W is 
negligible compared with AU. 

Overhauser has suggested that AU can be reduced by increasing m. 
While it is not convenient to analyse all possible combinations of 
spin-fluctuation states, it can be seen from (1) that increasing m 
reduces W exponentially, while AU is decreased only like m~*. Thus 
it appears that unless n is very near | either the non-magnetic or 
ferromagnetic Hartree-Fock states will have lower energy than 
Overhauser’s spin-fluctuation states. This conclusion is apparently 
valid for fairly general short-range repulsive interactions (including 
screened Coulomb interactions). 

Recently Kohn and Nettel (1960) have arrived at conclusions similar 
to ours in the weak interaction limit (n>1), by very general arguments 
which apply explicitly to the Hartree state. While our arguments are 
more particular, they show directly why exchange generates spin density 
fluctuations in the Hartree-Fock states for one but not for two or three 
dimensions. 
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The Etching of Sub-structures in Beryllium 


By J. Sawxitn and J. E. Mereprru 


Tube Investments Research Laboratories, Hinxton Hall, Cambridge 
[Received July 27, 1960] 


In the course of some work on the electropolishing of beryllium, using 
the commercial ‘ Disa Electropol’ apparatus and electrolytes, it was 
found that beryllium was being etched in addition to being electro- 
polished. With commercially pure cast beryllium both grain boundaries 
and sub-boundaries were delineated, and in addition there were etch 
figures in the grains. It was subsequently found that the commercial 
A2 electrolyte, containing, according to Jacquet (1956), perchloric acid, 
ethyl alcohol, and butyl cellosolve, could be used simply as an etchant 
on surfaces previously electropolished or mechanically polished. However, 
a fresh solution did not etch the boundaries. Only after electropolishing 
beryllium with it for a few hours did the solution work as an etchant. 
It was also discovered that the etchant did not pit the surface, but rather, 
it produced a deposit at the boundaries. For example, the boundaries 
could be ‘moved’ by brushing the surface, and observations in phase 
contrast illumination showed the boundaries to be standing proud of 
the general surface. 

An attempt was made to analyse the deposit by evaporating a carbon 
film on to an etched surface, stripping it, and examining the film, with 
the deposit attached, by electron diffraction in an electron microscope. 
The deposit could easily be recognized in the microscope, but it did not 
give a crystalline diffraction pattern. Recently, however, etched surfaces 
have been examined in an electron probe x-ray microanalyser and the 
deposit has been shown to be copper. 

The action of the etchant is thus to decorate the boundaries and other 
sites with copper, and it is presumed that the fresh solution did not etch 
because it did not contain copper ions, whereas the electropolishing of 
beryllium allowed the copper, present as an impurity in the metal, to 
go into solution. The next step was to add a few crystals of copper 
sulphate to the fresh A2 solution, and this then etched perfectly well. 
Finally, it has been shown that the A2 solution is not necessary as an 
ingredient. The best etchant found so far is an aqueous solution of 
copper sulphate, with the copper content at a few hundred parts per 
million. It is also possible to deposit silver at the boundaries using a 
silver nitrate solution, and gold using a gold chloride solution, although 
in the latter case the chloride ions cause excessive pitting of the surface. 
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Figures 1 and 2+ show respectively the microanalyser electron image 
and x-ray image taken in copper radiation from a surface etched in 
copper sulphate solution. Figure 2 shows conclusively the copper at 
the three grain boundaries and additionally at sites within the grains. 
It was possible to prove that the radiation came from the deposit and 
not, for example, from the grain boundaries themselves. 

Examples of the use of the copper sulphate etchant are given in figs. 3 
and 4. Figure 3 is an optical micrograph of cast beryllium showing 
sub-boundaries and grain boundaries. Figure 4 shows the complex 
polygonized structure produced by deforming hot extruded beryllium 
at room temperature and annealing for 24 hours at 850°c. In both cases 
the misorientations across the sub-structures were such that they could 
not be detected by polarized light methods. Some recent experiments 
on the heat treatment of beryllium indicate that impurities in the metal 
may play a part in the etching reaction. 
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Stress Concentrations at Crystal Surfaces and the 
Embrittlement of Sodium Chloride 


By D. M. Marsx 


Tube Investments Research Laboratory, Hinxton Hall, Saffron Walden, 
Hssex 


[Received July 4, 1960] 


THE brittleness of ionic crystals has attracted considerable attention 
recently and the purpose of this letter is to suggest an explanation of this 
effect in terms of stress concentrations at surface steps and other irregu- 
larities. The author has been studying the stress concentrations which 
may exist at the root of such surface steps using frozen-stress photo- 
elastic techniques on large models. The experimental results show that a 
step, as in fig. 1, exhibits a stress concentration factor at its root given by 


shear s.c.=tensile s.c.=1+0:74/(I/r). 


This type of relation has been verified over large range of values of J/r 
and for two different riser angles. The point of note is that the stress 
concentration is found to be almost equal to that at a Griffith crack of 
comparable dimensions, and so may be very severe. 


Fig. 1 


RADIUS r_ 


It is suggested that such stress concentrations embrittle ionic crystals by 
two alternative mechanisms: 


(1) In a ductile crystal, failure is initiated at an internal defect which is 
only effective as a crack source after severe work-hardening and deforma- 
tion of the whole crystal, as only then will it sustain a fatal stress 
concentration (Stokes etal. 1959a). However, a surface stress concentra- 
tion provides an alternative crack source which is available from the 
beginning of stressing and so may cause almost immediate failure. 
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(2) The surface stress concentration raises the stress locally so that a 
small length of crystal deforms and fractures before there is much overall 
strain. This embrittlement is assisted by the flow-induced fracture 
mechanism (Stokes et al. 1959 a) specific to ionic crystals, which limits the 
permissible local flow before fracture. . ; 

The experimental evidence for the first mechanism 1s provided by the 
fracture faces of many brittle specimens on which the surface defect 
responsible for failure can be clearly seen at the source of fracture. This 
experimental evidence was suggested by some interesting results of Stokes 
et al. (1959 b). They showed that NaCl single crystal bars could be made 
reliably and permanently ductile by washing, if every trace of solution 
was subsequently removed from the surface. However, their as-cleaved 
specimens were often brittle, and ductile specimens could be made brittle 
by two treatments :— 


(1) If a drop of solution was allowed to dry on the surface to form a 
‘ stain ’. 
(2) If carborundum was sprinkled onto the specimens. 


They further showed that while fracture of the ductile bars was initiated 
internally, the failure of brittle specimens started at the surface. However, 
the present author has found that not only is failure initiated at the surface 
in such specimens, but also that a surface step or notch could be found at. 
the source of the fracture. For example, fig. 2 (a) (Pl. 151) shows that 
fracture started from a step on the surface of a * stained’ specimen and 
fig. 2 (b) shows initiation at a cleavage step which intersects the fracture 
face. The effect of cleavage steps may be further demonstrated by 
deliberately cleaving specimens slightly skew to the {100} faces to give 
gross steps. Such specimens are invariably brittle. Fracture faces of 
carborundum-treated specimens show that failure starts from the roots 
of sharp notches which have been chipped out of the specimen surfaces. 
In contrast fig 2 (c) shows the internal initiation of fracture in a washed 
specimen. 

In general a complete range of surface-initiated fractures may be found 
from those in which fracture is almost completely brittle as in fig. 2 (b) to 
those in which the surface stress concentration is small and considerable 
flow has occurred before failure. In the latter instances it is presumed 
that the small surface stress concentrations have been assisted by the 
long-range stress concentrations, created by plastic flow, which cause the 
characteristic ‘ chequerboard ’ patterns observed in polarized light. 

Experimental evidence for the second embrittling mechanism (i.e. 
localization of strain) is provided by an experiment in which small 
semi-circular notches were washed out of NaCl bars before annealing. 
These were found to be brittle in four-point bending despite being given 
the standard washing treatment (Stokes et al. 1959 b) after the anneal. 
Micrographs in polarized light showed intense deformation at the notch 
only, while the fracture at the notch showed the characteristic ductile 
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failure as in fig. 2 (c). Thus premature failure was by localization of slip 
rather than by limitation of slip. 

One advantage of the present theory is that it provides a common 
_ mechanism for all of the above specimens and no evidence was found in 
NaCl for the special embrittling mechanism suggested by Stokes ef al. 
(1959 b) to account for the brittleness of their stained specimens. How- 
ever, microcracks such as they found in ‘as-cleaved’ MgO, certainly 
assist embrittlement when present. 

It is important to ensure that the proposed mechanisms are not equally 
applicable to metals as these do not exhibit this embrittlement. It is also 
important to ask why the surface steps and notches do not render them- 
selves harmless by flow as do stress concentrations in metals. Two 
obvious reasons are : 

(1) The surface stress concentrations can be more easily and quickly 
relieved in metals owing to the higher mobility of dislocations and the 
larger number of slip systems available. 

(2) In ionic crystals relief of the stress concentration by flow merely 
promotes failure by the second mechanism—a mechanism which is much 
more effective in ionic crystals due to the internal defect mentioned above. 

In summary, it has been suggested that the embrittlement of NaCl is 
often due to the presence of surface steps and notches which act as stress 
raisers and cause premature failure by two distinct mechanisms. A fuller 
and more rigorous treatment of these ideas will appear elsewhere. 
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REVIEWS OF BOOKS 


Elements of the Kinetic Theory of Gases. By E. A. GuaGENHEIM. (Pergamon 
Press, 1960.) [Pp. xii+92.] Price 17s. 6d. 


THE physicist or physical chemist whose knowledge of the kinetic theory of 
gases stops at about the level of the formula »=4nmél is often daunted by 
the mathematical complexity which accompanies more exact treatments 
of transport phenomena. Professor Guggenheim does not eschew mathematics 
—that is not his purpose—but by taking as his standard example the hard 
sphere model of a gaseous molecule he is able to explain the physical basis 
of the sophisticated theory with a minimum of computational involvement. 
Moreover, he shows that even this crude model can be made to give fairly 
close agreement with experiment when it is worked out correctly. The moral 
of this analysis, that this sort of data does not give any exact indication of 
the form of intermolecular forces, he drives home by a similar careful discussion 
of the second virial coefficient in terms of the scarcely less crude square well 
potential. In a final chapter he employs data for solid as well as gaseous 
argon to synthesize an interaction potential for this particular atom. Here, 
as in the rest of the book, the clarity of the presentation gives great pleasure ; 
one can recommend the book with enthusiasm to anyone who wishes to take 
his understanding of kinetic theory beyond the scope of a normal undergraduate 
course. ASSP? 


An Introduction to Thermonuclear Research. ALBERT Simon. (Pergamon 
Press, 1959.) [Pp. ix+182.] 35s. 


Tue full title of this book, which is An Introduction to Thermonuclear Research.— 
A series of lectures given in 1955, gives a very good indication both of its scope 
and limitations. Although an enormous number of papers on controlled 
thermonuclear research have appeared since the removal of security restrictions 
in 1958, there have been very few attempts to provide books which assemble 
the material in an orderly fashion for the benefit of the non-specialist reader. 
Dr. Simon’s book sets out to explain to such a reader what is being attempted, 
but does not attempt to explain what progress has been made in recent years. 
For no attempt has been made to refer to material published since 1956, and 
the reader who is looking for an up-to-date account of the subject will be 
disappointed. But the book does provide an excellent introductory account 
of the more important proposals for making thermonuclear reactors. It is 
clearly written, and the author explains the subject by means of simplified 
qualitative arguments, avoiding long and complicated mathematics. 
The early chapters deal with the potentialities of fusion reactors, the basic 
conditions which have to be met in such reactors, and the motion of charged 
particles in magnetic fields. In the next four chapters, the author explains 
the fundamentals of the Stellarator, of mirror machines and of pinch devices. 
Finally there is a chapter on stability and another on the diffusion of plasma 
across magnetic fields. The emphasis throughout is on the principles involved, 
and very little is said about experimental apparatus or results. Le Be 
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ERRATA 


On Dislocations Formed by the Collapse of Vacancy Discs, by C. ELBAUM, 
1960, Phil. Mag., 5, 669. 
Page 670, line 12: increases should read decreases. 


Page 671, eqn. (5): the equation should read as follows: 
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E. ROBERTS and R. W. K. HONEYCOMBE Phil. Mag. Ser. 8, Vol. 5, Pl. 146. 


Fig. 1 
STRESS 
STRAIN 
(a) (6) (c) 


Hysteresis loops from an aluminium crystal fatigued at +700 g/mm? resolved 
, shear stress at room temperature. (a) 400 cycles (0-6% life), (6) 2800 
eycles (4% life), (c) 40 300 cycles (57° life). 


teresis loops from a crystal fatigued at +1015 g/mm? resolved shear 
Hys saree a —78°c. (a) 6500 cycles, (6) 6500 cycles, rested at room 
temperature then re-tested, (c) 6600 cycles. 
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The M emission band of solid potassium. 


Kaos GRISP Phil. Mag. Ser. 8, Vol. 5, Pl. 148. 


Emission lines from solid potassium in the 40-1000 A range. 
Upper left: Ist O.L,,V, Ist O.LI, Ly, Ist OLX. 
Upper right: 2nd O.LI, Ly, 2nd O.X, Ist O.Y. 
Lower left: 4th O.Ly, 4th O.LI, 4th O.X. 
Lower right: Ist O.L,1. 
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Microanalyser electron image of junction of grain boundaries. 
Etched in copper sulphate solution. x 1450. 


1 . » . . r 
Same field as fig. 1, x-ray image taken in copper Ka radiation. x 1450. 
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Optical micrograph of cast beryllium showing grain and sub-boundaries. 
Etched in copper sulphate solution. x 55. 


Fig. 4 


Complex polygonized structure in deformed and annealed beryllium. 
Etched in copper sulphate solution. x 110. 
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(a) ‘Stained’ specimen (x 570). 


(b) As-cleaved specimen (x 100). 


(c) Standard washed specimen (x 100). 


Fracture faces of NaCl. 


